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FOREWORD 


This  report,  OSURF  Report  Number  2183-4,  was  prepared  by 
The  ElectroScience  Laboratory,  Department  of  Electrical  Engineering, 
The  Ohio  State  University  at  Columbus,  Ohio.  Research  was  conducted 
under  Contract  AF  19(628)-5929.  Dr.  John  K.  Schindler,  C-RDG,  of  the 
Air  Force  Cambridge  Research  Laboratories  at  Bedford,  Massachusetts 
was  the  Program  Monitor  for  this  research. 

The  material  contained  in  this  report  is  also  used  as  a  thesis 
submitted  to  the  Department  of  Electrical  Engineering,  The  Ohio  State 
University  as  partial  fulfillment  for  the  degree  Master  of  Science. 


ABSTRACT 


A  ray-fixed  coordinate  system  is  introduced  and  used  to  derive 
a  nsw,  compact  form  of  the  dyadic  diffraction  coefficient  for  an 
electromagnetic  wave  incident  on  a  perfectly-conducting  wedge. 

This  diffraction  coefficient  is  merely  the  sum  of  two  dyads; 
furthermore,  with  the  use  of  simple  correction  factors  which  have 
the  same  form  for  plane,  cylindrical,  conical  or  spherical  waves 
incident  on  the  edge,  the  dyadic  diffraction  coefficient  is  valid 
in  the  transition  regions  of  the  shadow  and  reflection  boundaries. 
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CHAPTER  I 
INTRODUCTION 

The  purpose  of  the  research  described  in  this  report  is  to  derive 
a  dyadic  diffraction  coefficient  for  electromagnetic  waves  incident 
on  a  perfectly-conducting  wedge.  It  is  desired  that  this  diffraction 
coefficient  be  valid,  not  only  in  the  illuminated  and  shadowed  regions 
surrounding  the  wedge,  but  also  in  the  transition  regions  adjacent  to 
the  shadow  and  reflection  boundaries.  Furthermore,  it  is  desired 
that  the  diffraction  coefficient  be  compact  and  accurate  so  that  it 
is  useful  for  computational  purposes.  These  objectives  have  been 
attained  in  the  form  of  a  new  dyadic  diffraction  coefficient,  which 
is  derived  in  the  pages  to  follow. 

The  diffraction  coefficient  is  introduced  in  the  geometrical 
theory  of  diffraction  in  order  to  relate  the  field  of  the  ray  in¬ 
cident  at  the  edge  to  the  field  of  edge-diffracted  ray.  The  dif¬ 
fraction  coefficient  is  analcgous  to  the  reflection  coefficient  used 
to  relate  incident  and  reflected  fields.  A  description  of  the 
geometrical  theory  of  diffraction  is  given  in  References  1,  2,  3. 

If  the  problem  involves  scalar  waves,  the  diffraction  coefficient 
is  a  scalar  quantity;  on  the  other  hand,  if  the  problem  involves 
vector  waves,  as  it  does  in  the  case  of  electromagnetic  waves, 
the  diffraction  coefficient  is  a  dyadic.  The  diffraction  coeffi¬ 
cient  in  question  is  found  from  an  asymptotic  solution  of  the 
pertinent  wedge  diffraction  problem,  where  the  result  can  be  given  a 
ray-optical  interpretation.  The  diffraction  by  a  wedge  has  received 
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a  great  deal  of  attention  over  the  years,  and  some  of  the  more  im¬ 
portant  contributions  to  its  asymptotic  solution  are  described  in 
the  following  paragraphs. 

MacDonald4'  obtained  an  integral  representation  for  the  eigen¬ 
function  expansions  of  the  total  field.  Plane,  cylindrical  and 
spherical  wave  illuminations  of  the  wedge  by  scalar  waves  were 
treated.  However,  the  earliest  asymptotic  solution  to  the  wedge 

r: 

problem  appears  to  have  been  due  to  Sommerfeld.  He  considered 
the  case  of  scalar  plane  wave  illumination,  and  his  solution  is 
valid  outside  of  the  transition  regions  at  the  shadow  and  re¬ 
flection  boundaries;  however,  it  is  not  valid  within  the  trarisi- 

fi  7 

tion  regions.  Pauli  and  Oberhettinger  obtained  uniform  asymp¬ 
totic  solutions  which  can  be  used  in  the  transition  regions.  Both 
of  these  solutions  are  restricted  to  the  exterior  wedge  problem, 
i.e.,  they  are  restricted  to  wedge  angles  less  than  180  degrees; 
furthermore,  there  is  a  deficiency  in  Pauli's  solution  that  im¬ 
poses  a  further  limitation.  This  deficiency  is  pointed  out  in 

8 

Reference  9.  Felsen  has  obtained  uniform  asymptotic  solutions 
for  plane,  cylindrical  and  spherical  wave  illuminations  of  the 
wedge.  His  solutions  can  be  applied  to  both  the  interior  and  ex¬ 
terior  wedge  problems.  In  addition  to  wedges  with  the  usual  Dirichlet 
and  Neumann  Boundary  conditions,  Felsen  has  considered  wedges  with 
impedance  boundary  conditions.  Pauli  obtained  a  factor  type  cor¬ 
rection  term  which  makes  his  solution  valid  in  the  transition  region, 
whereas  Oberhettinger  and  Felsen  have  additive  type  corrections. 

These  differences  result  from  the  different  methods  used  to  treat 
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the  pole  close  to  the  saddle  point  in  the  asymptotic  solutions. 

9  14 

Hutchins  and  Kouyoumjian  *  generalized  Pauli's  solution  for  the 
scalar  plane  wave  illumination  of  the  wedge  and  obtained  a  factor 
type  correction  for  the  transition  regions,  which  makes  the  leading 
term  in  their  asymptotic  expansion  more  accurate  than  those  of  the 
solutions  described  previously.  Their  solution  corrected  the  afore¬ 
mentioned  deficiency  in  Pauli's  solution.  In  addition,  the  relatively 
compact  form  of  their  solution  further  adds  to  its  desirability  from 
the  computational  point  of  view. 

The  work  described  thus  far  has  been  limited  to  the  case  of 

scalar  waves.  The  more  difficult  vector  diffraction  problem  has 

been  treated  by  Nomura,*0  who  obtained  a  uniform  asymptotic  solution 

for  the  problem  of  an  oscillating  electric  dipole  illuminating  a 

perfectly-conducting  wedge.  Even  though  his  solution  employs  a 

factor  type  correction  for  the  transition  region,  the  expressions 

for  the  field  possess  a  rather  complicated  form  and  they  are  not 

11  12 

amenable  to  easy  computation.  Tuzhilin  *  treated  the  scalar 
and  vector  wedge  diffraction  problems;  his  results  possess  additive 
type  correction  terms  for  the  transition  region  which  involve  com¬ 
plicated  integrals  for  the  cylindrical  and  spherical  wave  illumi- 

13 

nations.  Recently  Bowman  and  Senior  unified  and  simplified  the 
results  of  Tuzhilin  and  others  by  treating  the  case  of  a  conducting 
half-plane  excited  by  an  oscillacing  electric  dipole.  An  asymptotic 
representation  for  the  Hertz  potential  was  presented  with  an  ad¬ 
ditive  type  correction  for  the  transition  region.  The  results  of 
Reference  13,  although  simpler  than  those  ot  References  11  and  12  and 
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others,  still  lead  to  cumbersome  expressions  for  the  fields  obtained 
from  the  space  derivatives  of  the  Hertz  potential.  In  summary,  none 
of  the  electromagnetic  solutions  described  in  References  10  through 
13  appear  to  be  in  a  form  suitable  for  the  derivation  of  a  compact 
dyadic  diffraction  coefficient. 

It  was  decided  to  derive  the  dyadic  diffraction  coefficient  for 
several  types  of  edge  illumination.  Although  the  diffraction  coeffi¬ 
cient  is  independent  of  the  type  of  edge  illumination  outside  of  the 
transition  regions,  it  differs  with  the  type  of  edge  illumination 
within  the  transition  region.  Let  the  edge  of  the  wedge  lie  along 
the  z-axis.  In  the  case  of  plane  and  cylindrical  waves  incident  on 
the  perfectly-conducting  wedge,  the  field  components  are  readily 
found  from  Ez  and  H7,  the  z-components  of  the  electric  and  mag¬ 
netic  fields,  respectively.  In  turn,  the  solutions  for  and  Hz 
reduce  to  the  ordinary  scalar  diffraction  problems.  The  incident 
spherical  wave  may  have  z-directed  electric  or  magnetic  current 
dipoles  as  its  source.  The  electric  and  magnetic  fields,  due  to 
these  dipoles,  can  be  obtained  from  F  and  A  ,  the  z-components  of 
electric  and  magnetic  vector  potentials,  respectively.  Again,  the 
solutions  for  Az  and  Fz  reduce  to  ordinary  scalar  diffraction 
problems. 

In  Chapter  il  asymptotic  solutions  of  the  scalar  diffraction 
problems  are  obtained  by  the  modified  Pauli -Clemmow  method  of 
steepest  descent.  The  wedge  illuminated  by  a  line  source  is  treated 
first  and  an  integral  representation  of  the  field  is  obtained  from 
its  eigenfunction  expansion.  By  letting  the  line  source  recede  to 
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infinity  this  integral  representation  is  transformed  into  an  integral 
representation  for  a  plane  wave  incident  on  the  wedge.  The  integral 
representation  for  the  field  of  the  point  source  illumination  is  de¬ 
duced  from  the  integral  representation  for  the  field  of  the  line 
source  illumination  by  a  Fourier  transformation  involving  the  z- 
coordinate.  After  asymptotically  evaluating  these  integrals,  the 
scalar  diffraction  coefficients  and  their  correction  factors  for 
the  transition  regions  are  obtained.  These  diffraction  coefficients 
are  of  interest  not  only  in  the  case  of  scalar  diffraction,  but  also 
in  the  case  of  vector  diffraction,  where  they  appear  in  the  dyadic 
diffraction  coefficient. 

In  Chapter  III  the  diffraction  of  electromagnetic  waves  by  a 
perfectly-conducting  wedge  is  treated.  The  plane  and  cylindrical 
wave  solutions  of  Chapter  II  are  generalized  to  oblique  incidence, 
so  that  in  the  place  of  an  incident  cylindrical  wave  one  has  an 
incident  conical  wave  whose  source  is  a  traveling  wave  electric  or 
magnetic  current  line  source.  Starting  with  Ez  and  Hz  one  then 
determines  integral  representations  for  all  of  the  field  components. 
These  integrals  are  then  evaluated  asymptotically. 

In  order  to  simplify  the  expression  for  the  diffracted  electro¬ 
magnetic  field,  special  coordinate  systems  for  the  incident  and 
diffracted  fields  are  introduced.  These  coordinate  systems  are 
suggested  by  the  law  of  edge  diffraction  and  they  are  fixed  in  the 
incident  and  diffracted  rays.  'When  the  components  of  the  incident 
and  diffracted  fields  are  written  in  this  coordinate  system,  and 
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only  the  leading  term  in  the  asymptotic  solution  is  retained,  the 
resulting  expression  for  the  dyadic  diffraction  coefficient  reduces 
to  a  sum  of  two  dyads.  If,  instead,  a  simple  edge-fixed  coordinate 
system  were  used,  the  resulting  dyadic  diffraction  coefficient  would 
appear  as  a  sum  of  seven  dyads*,  thus,  the  importance  of  introducing 
the  special  ray-fixed  coordinates  can  hardly  be  overemphasized. 

The  case  of  a  spherical  wave  obliquely  incident  on  a  perfectly- 
conducting  wedge  is  treated  in  a  similar  manner;  in  this  case,  the 
field  components  are  found  from  Az  and  F  . 


CHAPTER  II 

SCALAR  WAVE  DIFFRACTION 


This  chapter  deals  with  the  scalar  v/ave  diffraction  by  an  im¬ 
penetrable  wedge;  both  Dirichlet  and  Neumann  boundary  conditions 
are  treated.  Plane  wave,  cylindrical  wave  and  spherical  wave 
illuminations  of  the  edge  are  considered.  The  relationship  of 
these  problems  to  the  electromagnetic  problem  is  indicated. 

A.  Cylindrical  Wave  Illumination 

Consider  a  z-directed,  uniform  line  source  of  unit  strength 
radiating  cylindrical  scalar  waves  in  the  presence  of  an  impenetrable, 
infinite  wedge  as  shown  in  Fig.  1.  The  total  field,  which  consists  of 


Fig.  1.  The  2-D  wedge  and  line  source  in  the 

circular  cylindrical  coordinate  system. 
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the  incident  and  scattered  components  surrounding  the  wedge ,is  given 

by  a  two-dimensional  scalar  Green's  function  G  ( p~,p*' )  where 

h 

(1)  (v2  +  k2)  GS(F»P')  =  »6(|r  -  ?!) 

h 

2 

=  the  tWG-dimensional  Laplaci&n  operator  in  circular  cylindrical 
coordinates,  and  6{  |pT  -  p"'  J )  is  the  Dirac  Delta  function,  k  =  wave 
number  of  the  linear,  homogeneous  and  isotropic  medium  surrounding  the 
wedge.  )  satisfies  the  following  boundary  conditions  at  the 

wedge  surface: 

(2)  Gs(F,p')  =  0  at  <}>  =  0,  rnr; 

(3)  —■  ,p * )  =  0  at  <j>  =  0,  mi. 

Boundary  condition  (2)  corresponds  to  the  acoustic  case  for  a  soft 
wedge  whereas  (3)  corresponds  to  the  acoustic  case  for  a  hard  wedge. 

In  mathematical  terminology,  (2)  corresponds  to  the  homogeneous 
Dirichlet  boundary  condition  and  (3)  corresponds  to  the  homogeneous 
Neumann  boundary  condition. 

Furthermore,  G  (p",pT' )  also  satisfies  the  Sommerfeld  radiation 

15  h  ig 

condition  and  the  Meixner  edge  condition,  respectively.  Thus 

the  solution  for  G{.(p’,p'‘)  is  unique, 
ft 

The  above  problem,  depicted  in  Fig.  1,  represents  the  basic 
canonical  problem  in  this  report.  It  will  be  used  to  deduce  the 
solutions  of  the  other  canonical  problems. 
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The  Green's  function  Gs(p",p 1 )  for  this  canonical  problem  is 

n 

given  by  a  convergent  eigenfunction  expansion  as 


m=0 


(4)  Gs(p,D  =  =  $  l  Ojkp)  h‘%p') 

h  h  1,1=0  K  n 


where 


[cos  ^  U-d>')  +  cos  ($+$')], 


0  <_  4>,  <j> '  <_  mr 

0  <  p  <  p'  <  ®  , 


and 


e  =  1,  m  =  0 
m  ? 

=  2,  m  ?  0. 

0  (kp)  and  H^(kp')  represent  the  cylindrical  Bessel  function  of 

F  F 

the  first  kind  and  the  cylindrical  Kankel  function  of  the  second 

kind,  respectively.  A  time  dependence  of  the  e’fJut  type  is 

assumed  and  suppressed.  represents  the  modal  description  of 

the  waves  existing  in  the  presence  of  the  wedge  structure  excited 

by  a  line  source  at  {p * ) •  £>s  ss  given  by  (4)  converges  for  all 

Pj'PjP1  »$'  and  n  (interchanging  p  and  p*  in  (4)  yields  Gs  for  p  >  p'), 

h 

It  is  easy  to  see  that  the  total  electric  field  (z»airected) 
is  proportional  to  Gg  if  the  line  source  is  electric,  whereas 
the  total  magnetic  field  (z-directed)  is  proportional  to  Gh 
if  the  line  source  is  magnetic.  Thus,  the  scalar  problem  is 
useful  in  the  treatment  of  the  electromagnetic  problem  of  a  uniform 
electric  or  magnetic  current  line  source  illuminating  a  wedge  which 
is  perfectly  conducting. 
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v 


If  the  line  source  is  elactric,  then  the  total  electric  field,  Ez 
is  given  by 

(5)  EzGT)  =  -jwjj  IGS, 

where  tu  is  the  angular  frequency  and  y  is  the  permeability  of  the 
medium  surrounding  the  wedge. 

If  the  line  source  is  magnetic,  then  the  total  magnetic  field, 
Hz  is  given  by 

(6)  HZ(F)  =  -jus  MG^, 


where  e  is  the  permittivity  of  the  medium  surrounding  the  wedge  and 
I  and  M  are  the  strengths  of  the  electric  and  magnetic  current  line 
sources  respectively. 

An  integral  representation  for  the  product,  dm(kp)H^2^(kp')  is^ 

rT  rf 


(7) 


Om(kp)H<2)(kp' 


)  = 


i- 


[c“j“  Jit-k^p^p'^t"1] 


n 


n 


c  >  0,  -jE  >  -1  and  |pj  <  jp'  |  , 


.2  . 


Im  (---  |P—  )  represents  the  modified  cylindrical  Bessel  function 


n  2 

of  the  first  kind.  A  useful  integral  representation  for  1^  ) 

.  4 

is  given  by 


m 

n 


(8)  I.  (^4  - 


2tt 


y'+J-  cos  C  +j  J  C 


dC» 


Y+j° 
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where  -it  <  y*  <  0, 
it  <  y  <  2  fr. 


where  the  contours  L  and  L'  are  indicated  in  Fig.  2.  Replacing  the 
cosine  terms  in  (4)  by  exponentials,  and  utilizing  (7)  and  (9)  in 
(4),  interchanging  the  orders  of  summation  and  integrations,  and 

1.2  i 

choosing  the  L'  contour  representation  for  l0(t~~)  (corresponding 
to  ib  =  0),  allows  one  to  write  (4)  as 


(10)  G$(p,pe)  =  G(p,p';e")  +  G(p’,p‘;3+), 

where 


and 

(12) 


B+  5  $  +  <J>'  =  P 
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Fig.  2.  The  £  and  L'  contours  for  the  integral  representation  of 
lm(- in  the  complex  5  plane, 
n" 


Eq.  (11)  may  be  further  simplified  by  noting  that 


*  j  jjf  (e+s)  “  i  ?  U+e) 

I  e  =  S  e  "  -!• 

m=l  m=0 


and 


(13)i 


"  ±3i(5+B> 

1  6 

m=0 


1 


1  -  e 


Also 


*  0  ~  U+3) 


m=0 


_  1  -  _  1_  rnt  (  111.  \  _ 

*  ^  '■  2n  ' 


2j 


1 

2  * 


and 


09  -0  "  (C+3)  1  c+o  ,  1 

J0  e  '  ?J  cot  (  sr  > +  2 
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Therefore,  (11)  may  be  written  via  (13)  as 


(14)  G(p,0-;e)  -  ■■  -f  fC'J"  .WWftpVV1]  . 

8’ »  l 

[Jj< ~(S*)-i]« 

which  on  interchanging  the  orders  of  integration,  gives 


5  + 


[.  dt 
t  * 


(15)  G(p,p';b) 


8ir2n  \ 


d? 

L  l 


2j  [  2n  j 


ll- 


L. 


fC-j“  dt  %[t-k2(o2+p'2)t”1l  + cos  i 


lx  {,  d5'  « 


fC-J° 


dt 


r  • 


>sLt-k2 ( p2+p  * 2)  t“ ■*■]  +  cos  5 


u0 


It  can  be  shown  thatH 


(16) 


“  2  J 


°  £  K0(  jz) » 


C-Jc 


where  KQ( jz)  is  the  modified  cylindrical  Bessel  function  of  the 
second  kind,  of  order  zero  and  argument,  jz.  Utilizing  (16)  in 
(15)  gives 
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or,  on  combining  terms  one  obtains 
(17) 


g(p,p‘ ;3) 


1 

4TrZn 


j  KQ(jk  v/p2+p'“-2pp' 
L-L' 


cose) 


Let 
then , 


k  /p2  +  p'2  2pp'  cos  c  “  zU) 


~V"  f  iKo(jz(5))d5  f  |K0(32(5))de  -  f  |  K0(jz(e:))c'c 

4ti  n  LlL,  [  |s 


8rr2n 


r 

' 

k  (jz(e))de  - 

K0(jz(-e))d(-e)  ' 

J 

U  -i 

-  </ 

2 

4tt  n 


f  =  0, 


since  z(e)  =  z(-e). 

Finally,  (17)  can  be  written  as 


(18a)  G(p,pl ;3+)  = 

4tt  n 


d<H 


L-L* 


K0(ok^+P--2PP-  cose)  « 


|2_ 


Equation  (18a)  directly  leads  to  a  complex  contour  integral  repre¬ 
sentation  for  (p,p')  =  G(p,p';e")  +  G(p,p';6+).  The  integral 
representation  is  still  an  exact  form  of  the  solution. 
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If  |zU)|  is  large  with  respect  to  the  order  of  K0(jz(s)), 
one  may  use  the  large  argument  approximation 

K0(jz(c))  *  J  237(77  where  z{?)  =  k^p+p,2-2pp‘  cos  5 

and  -n  arg  ( jz)  <_  v. 

In  the  asymptotic  solution  described  later  it  will  be  seen  that  |z| 
remains  large  in  the  neighborhood  of  the  saddle  points.  Therefore, 
(18a)  may  be  written  as 


(18b)  G(p,p’;f3T)  f  d5  /— 

Stj  jn  LlL,  J2j 


2jk(p^4-p  1  ^-2pp  *  cos  K)h 
f  \  «jk/ p2V2-2pp'  cos  ? 


•  cot  (fw6-)  8 

Furthermore,  the  exponential  may  be  approximated  as 


(19)  e 


-jk/ p2+p,2-2pp'cos5  _  0-jk/  (p+p')2-2pp'(cos  5  +  1) 


-jk(pV)  U 

ts,  -e 


- sr  (  1  +  COS  ?)  } 

(p+p') 


The  above  approximation  will  be  justified  a  little  later  in  the 
chapter. 

Let 


FX(§,P+)  = 


_i_  / . . l _ 

8Tr2jn  J  2jk(o2+p'2-Rpp'  cos  fJ% 


e-jk(p+p') 


and 


f(c)  =  j  [1  +  cos  Cl  . 
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Thus, 

(20)  6(p.p- ;e*)  *  f  d5  r^J)  e“f(5) 

L-L' 

is  a  compact  notation  for  the  expression  in  (18b)  but  containing  the 

approximation  of  ( 19). k  =  and  is  assumed  to  be  large.  Thus, 

Eq.  (20)  is  in  the  proper  form  for  an  asymptotic  evaluation  of  the 
integral  via  the  method  of  steepest  descent.  The  steepest  descent 
paths  (SDP)  selectfc  for  the  asymptotic  evaluation  of  (20)  are  those 
which  pass  through  the  saddle  points  of  f(c)  at  ?s  =  ±?r  for  the 
following  reason. 

The  saddle  points  of  f(§)  occur  at 
af(  , 


0; 


however  5,  -  ±*are  the  only  ones  chosen  because  the  steepest  descent 
paths  through  £s  =  ±*  allow  one  to  close  the  (L- L‘)  contour,  thereby 
facilitating  the  use  of  the  Cauchy  Residue  Theorem.  Figure  3  indicates 
the  locations  of  the  steepest  descent  paths  through  the  pertinent  saddle 
points  at.  =  ±n.  It  is  easily  seen  that  the  required  integral  in 
(20)  ’.s  given  by: 


(21) 


d£  FU,3+)  e*f^d£  =  -  j  d£  F(5,b?)  eKf^5'  dg  - 
L-L'  SDP(ir) 


3*)  e^U)  + 


j  dsFjU* 

SDP(-ti) 

sum  of  the  residues  of  FtU,.3+)  enclosed!  + 
.between  L-L1  and  SDP(±tt)  -1 


+  2uj 


+  (branch  cut  contributions,  if  any)  . 
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In  order  to  ensure  convergence  of  the  integral  in  (18b),  one  must 
have  Re  f(c)  <  0  as  jlm  £|  -*•  ».  The  portions  of  the  complex  £  plane 
for  which  Re  f(0  <  0  as  |Im  £|  ■+  »  are  shown  as  shaded  vertical 
strips  in  Fig.  3.  The  SOP  equations  are  found  from  the  condition 
that 

Int  fU)  =  Im  f(£s) 
or 

cos  u  cosh  v  =  -1,  where  £  =  u  +  jv. 


The  above  equation  is  valid  for  £g  =  ±ir  and  the  paths  chosen  are 
depicted  in  Fig.  3,  the  choice,  as  mentioned  before,  being  governed 
by  the  fact  that  Re  f($)  <  0  on  the  SDP(±rr)  as  |Im  £|  -»-  »  (where 
Im  £  =  v),  and  that  these  paths  blend  in  with  the  L  -  L'  contours  as 
| Im  £ |  -*-«>,  thereby  allowing  one  to  obtain  a  closed  path  of  in¬ 
tegration.  (The  paths  corresponding  to  the  images  of  SBP(±n)  about 
Re  £  =  ±n  also  satisfy  cos  u  cosh  v  =  -1  but  are  paths  for  which 
Re  f(£)  >  0  as  |  Im  £ |  -*-«»»). 

It  is  desirable  at  this  juncture,  to  study  the  singularities 
of  FU.e*)  in  the  complex  £  plane. 


r(£,B+)  = 
1 


i  /  .  * - rCOt(^Vjk(p+p,). 

8u2jn  J  2jk(p2+p,2-2pp'  cos  £)2  \  n  / 


F(£,a:t)  has  branch  point  singularities  at  £  =  £b  which  are  the  roots 
1 

of 

P2  +  p'2  -  2pp 1  cos  £b  =  0 
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For  the  case  where  p  f  p',  it  follows  that  p  +  p'  >  2pp* ,  and 
therefore 

-1  2  +  '2 

=  ?.% n  ±j  cosh_i  p-  2pp,f? —  *  where  i  =  0,  ±1,  ±2,  ••• 

The  branch  points  and  the  branch  cut  locations  are  indicated  in 
Fig.  3,  where  it  is  seen  that  the  branch  points  are  not  in  the 
vicinity  of  the  saddle  points  at  ?s  =  ±tt,  and  the  branch  cuts  lie 
outside  the  closed  path  formed  by  /.  -  L'  and  SDP(±ir),  so  that  the 
closed  path  is  unaffected  by  the  branch  cuts.  The  integrand  is 
single  valued  on  and  within  the  closed  path,  and  since  no  deformations 
of  this  path  are  necessary  for  single  valuedness,  the  branch  cut 
does  not  contribute  to  the  integral. 

The  pole  type  singularities  of  F,U)  occur  at  5  15  5n»  where  5 
are  specified  by 

Sp  =  -<3+  +  2nN*,  N  =  0,  ±1,  ±2,  ••• 

The  residues  corresponding  to  are  evaluated  for  |5p|  <  1  as  only 
these  poles  lie  within  the  closed  path  of  integration  since  e+  = 

$  +  <j>',  (it  is  easily  seen  that  Sp  is  a  real  quantity).  It  can 
be  shown  that  sp  are  simple  poles. 

For  sufficiently  large  tc,  the  major  contribution  to  the  in¬ 
tegrals  evaluated  over  SDP  (±tt)  occurs  only  from  the  "immediate 
vicinity"  of  the  saddle  points  (c$  =  ±tt);  hence  it  follows  that  the 
approximation  of  (19)  is  justified,  as  j(l  +  cos  ?)|  is  small  for 

The  first  order  saddle  point  approximation  for  a  pole  not  close 

* 

to  the  saddle  point  is  given  by 
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(22) 


SDP  (  ±tt) 


'  5U.6VfU)d5  -  r(±„,e*)  eKf(±',)  j 


u  u/4> 

,3^37r/4i 


i'r)  .. 


J  x  JV  4  > 

in  which  e  is  associated  with  j;s  =  7^ and  e  with  =  -n. 

The  residue  contribution  to  the  integral  over  the  closed  path, 

c  is  given  by 


(j)  FU,e+)fcKf^d5  =6  dg  =  2ttj  l 

7  1  7  QU.e+) 


pUg^£) 


Q'(5p.6T) 


where  c  =  (L  -  L')  +  SDP(tt)  +  SDP(-n),  and  P(s,P+)  and  QU,B+)  are 
both  analytic  functions  of  the  complex  variable,  5. 

-ik/2V2-Z 


P(5,3+)  = 


2>ik(p2+p'2”2pp'  cos  c)J 


cos(ii-) 


-jk/p  V  -2pP'cos  5 


and 


with 


QU,8+)  =  8iT2jn  sin(^~  j 
Q'U.S*)  =  4*2j  cos(^)  . 


The  residue  at  e^  =  -  B+  +  2nNn  is  thus  given  by 


(23) 


2irj 


.  I 


IL 


itk(p2+p '  2-2pp  1  cos  [-p++2nNu] )  ** 


.«,/ 2*.  .2. 


_ _  .  _  __9 _ *  ^ac  r-o  4-9nW— 1  !•  ,  4- 

jiNffj  *-pp  . 


z\ 


M  . . .  j  U[ir-|-p-+2nNir|] , 


which  may  be  recognized  as  the  asymptotic  form  of 

2^J  i  Ho2) (k^2+P ' 2"2pp ' COS [-e++2nN7T] )  '  U[ir-|-8++2nN*|] » 
and  therefore  represents  the  geometrical  optics  contribution  in  terms 
of  the  line  source  fields  due  to  the  source  or  it's  image  (the  total 
geometrical  optics  contribution  being  a  superposition  of  the  incident 
and  the  reflected  fields). 


U(t)  = 


0}  if  t <  0 
hi  if  t  =  0 
{  1,  if  t  >  0, 


U(0)  =  h  ensures  the  Cauchy  principal  value  when  the  pole  singularity 
lies  on  the  closed  contour  c.  The  function  U(Tr-|-3*+2nNTf|)  auto¬ 
matically  restricts  the  geometrical  optics  contributions  to  aspects 
defined  by  |-e*  +  2nNirj  £  ir,  so  that  the  contributions  from  poles 
outside  the  closed  contour,  c,  are  not  included. 

Thus,  for  the  N  =■  0  case,  one  obtains  for  the  geometrical  optics 

contribution,  G-?*0,  the  following  field 
h 


(24)  G?*°* 


.  1 
4 


ii 


nk(p  +p'  -2pp' cos (<{>-c> '))^ 


•  e-jk(pV2-2pp'cosU-*')),s  U(H„,|)}  ; 


t 


i  /  2  i  _ 

irk(p ^+p  1  ^-2p p  1  cos (' $+$ '  yfi 

.  6-jk(p£‘+p|2-2pp'cos(^+4i'))J2 


uU-U+*,|)f  . 
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(The  minus  (plus)  sign  between  the  two  terms  on  RHS  of  (24)  cor¬ 
responds  to  the  soft  (hard)  boundary  condition).  The  first  term  on 
the  right  of  (24)  is  the  direct  source  (at  p '  »<f> * )  radiation  at  the 
observation  point  (at  p,<}>),  whereas  the  second  term  represents  the 
reflected  wave  contribution  at  the  observation  point  (at  p,^)  which 
appears  to  emanate  from  a  virtual  line  source  (at  p  1  1 ) .  The  re¬ 

flection  occurs  from  the  wedge  face,  =  0.  In  general,  the  incident 
field  is  described  by  N  =  0  and  g"  =  4>-<t> * .  The  reflected  field  from 
the  surface  <p  =  0 resulting  from  the  incident  field,  is  described 
by  N  =  0  and  g+  =  $  +  Other  values  of  N  describe  fields  which 
may  be  reflected  either  from  the  surface  <j>  =  0  or  the  surfaces  <|>  =  nir. 

Since  the  total  field  is  composed  of  the  geometrical  optics  field 
and  the  diffracted  field,  it  follows  that  the  saddle  point  contri¬ 
butions  must  yield  the  diffracted  ray  field  in  accordance  with 
Keller's  theory.  The  first  order  saddle  point  results  for  large 
(kpp')/(p+p')  as  obtained  from  (22)  is  given  as 


Gs  =  ~ 


or, 


3DP(tr) 

r 

SDP(tt) 


F1U,e")e,cf(5,e  +  |  F1(5,6~)eKf(c,e  ^  + 


SDP(-n) 


FjU^e*1^’3  *d$  +  |  F1(5,B+)eKf^’8  V 


6«  .  £^1  j 

^  87rnk/pp'1 


r 

V 

e 


J2 


SDP(-tt) 

•  TT 


jf 

e  ^  cot 


cot(£iLj  +  c-J2cot(^.jj 
(#)+e'if  «*(&) 
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f.e. 


(25) 


d  i  e"J^(p+p' ) 
GS  ^  ‘SfjK  r— p 

h  ‘'pp 


1  IT 

Wsln  n 


cos  f  -  cos  ■— 
n  n 


1  C1-  tr 

~  sm  “ 
n  n 


cos  ~  -  cos  — 
n  n 


For  N  =  0  case,  the  total  field  surrounding  the  wedge  is  given  asymp¬ 
totically  by  superposing  the  results  of  (24)  and  (25)  as  indicated 
by  (21). 

The  value  of  the  incident  wave  at  the  edge  is  given  by,  u1  where 

u’  “  -jj-H<2W)  '  -  i  iH-jr  c'jkp,> 

for  large  kp'.  Therefore,  in  accordance  with  the  geometrical  theory 


of  diffraction,  may  be  written  as 
h 

(26)  Gs  '  [u1  '  - 


2“jkp 

~ 


for  large  kpp,/(p+pl),  where  Ds  is  the  wedge  diffraction  coefficient 
associated  with  either  the  soft  or  the  hard  boundary  condition  at 


the  wedge  surface.  The  factor 


,-jkp 

/T 


describes  the  variation  of 


the  amplitude  and  phase  of  the  field  of  the  diffracted  rays  emanating 
from  the  edge  and  traveling  in  the  positive  p  direction. 

Comparing  (25)  with  (26)  it  is  seen  that 


(27) 


_  e-JTT/4 

1  TT 

—  sin  — 
n  n 

-cos  n  - 

mu  i  ta  lcu  i 

-  /  mX  j  j _ _ i 

ii  \ct )  id  lucnuua  i 

irsi"K 


cos  ^  -  cos 


Sommerfeld.  This  result  fails  at  the  shadow  and  reflection  bound¬ 
aries  where  +  2Nu,  and  it  is  evident  that  at  the  shadow  and 
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reflection  boundaries  the  diffraction  coefficient,  becomes  singular. 

As  mentioned  earlier,  <j>  =  4> 1  +  ir  (N=0)  defines  the  shadow 
boundary;  however  a  shadow  boundary  does  not  occur  if  the  exterior 
wedge  angle  nu  <  ir  or  if  n-rr  ~  7r  <4>*  <  n  in  the  case  where  nn  -  tt>  0. 

A  reflection  boundary  associated  with  the  wedge  surface  4>  =  0  occurs 
at  <p  =  tt  -  4> *  (N=0);  if  mr  >  tt  -  <j>'.  Other  reflection  boundaries 
may  result  from  the  reflections  at  the  two  surfaces. 

In  the  transition  regions  of  the  shadow  or  a  reflection  boundary 
a  pole  is  close  to  the  saddle  point.  Later,  we  will  be  interested 
in  the  pole  closest  to  the  v  saddle  point,  and  the  pcle  closest  to 
the  -n  saddle  point.  The  values  of  N  associated  with  these  poles 
are  defined  to  be  N+  and  N“,  respectively.  The  value  of  N*  is 
decided  by  the  integer  which  most  nearly  satisfies  the  equation: 

2nN±Tr  -  6  =  ±tt  , 
where  e  =  3+  =  4>  +  ^ 1 . 


Whenever  the  poles  of  the  geometrical  optics  fields  lie  in  the 
vicinity  of  the  saddle  points  which  happens  when  the  observation 
point  is  in  the  transition  region,  the  ordinary  steepest  descent 
method  is  not  applicable  and  a  more  sophisticated  approach  is 
necessary.  This  approach  accounts  for  the  pole  singularity  of 
F^U)  when  it  is  in  the  vicinity  of  the  saddle  points  of  fU);  and 
this  method  will  be  referred  to  as  the  "Pauli-Clemmow  modification  of 
the  rnethod  of  steepest  descent11.  *  *  mis  wouii  ieu  steepesc  descent. 


method  has  been  presented  for  a  general  case  with  sufficient  details  in 
Appendix  I.  Hence,  no  formal  derivations  will  be  presented  for 
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the  different  problems  corresponding  to  the  different  source  il¬ 
luminations  considered  in  this  report.  Only  the  results  will  be 
presented  for  the  sake  of  brevity.  A  summary  of  the  details  in 
Appendix  I  may  be  stated  in  the  following  manner. 

The  modified  steepest  descent  method  essentially  involves 
replacing  the  integrand  by  a  product  of  two  functions,  one  of 
which  is  singular  with  a  simple  pole  and  describes  the  effect  of 
the  pole  near  the  saddle  point  in  question;  the  other  factor  is 
analytic  in  the  neighborhood  of  the  pole  and  the  saddle  point 
in  question.  The  analytic  portion  of  the  integrand  is  then  ex¬ 
panded  in  a  MacLaurin  series  about  the  saddle  point  and  the 
resulting  series  is  integrated  termwise.  In  all  of  the  results 

r 

presented  in  this  report,  only  the  first  term  of  the  MacLaurin  ex¬ 
pansion  is  retained  for  simplicity.  Thus,  the  results  presented 
here  correspond  to  a  first  order  asymptotic  approximation.  Equation 
(25)  is  corrected  through  the  modified  steepest  descent  method  of 

Appendix  Ijfrom  where  (diffracted  fields)  for  large  k  =  (kpp')/(p+p‘) 

h 

is  directly  given  by 


(28)  Gg  * 

h 


a-jk(p+p') 

1  *  1  2ti_(p+p ’.)  [| 

r  ii 

/ 

,e  cot  ^ 

'«+b"\ 

rsr) 

F[ica+(3  )]" 

8n2jn  ' 

j2jk(p+p'F  J  <ppr  y 

.3tt 

e  4  cot  ^ 

F[<a-(ff)]|  + 

i11 

Jt 

■  e  cot 

( ft+&+' 

{  “ST. 

)  F[Ka+(e+)]- 

,.3tt 

J 

> 

e  ^  cotj 

F[<a  (6  )]jj  . 
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b 

c 

t 

l 

t 


^  j 

(Note:  g  =  4>  +  <P').  As  before,  G<?  may  be  re-expressed  for  inter- 
pretation  in  terms  of  the  GTD  as  follows. 


(29) 


's 

h 


J 


where,  u1  represents  the  line  source  field  incident  at  the  edge.  The 

new  diffraction  coefficient,  denoted  by  D  is  given  as 

h 

(30)  Ds  =  (d+(g~,n)  F[>a+(g_)]  +  d'(e",n)  F  Ka“(e")]>  + 
h 

+  (d+(6+,n)  FCKa+(g+)]  +  d“(g+,n)  FCKa"(e+)]}  , 

where 

<31)  r(e,n)  5 '  tis  * cot  (^)  • 


2  =  <  =  kpp!/(p+p') 

and 


where  the  positive  branch  of  the  square  root  is  taken,  and 


(33)  a^g)  =  (1  +  cos(-g  +  20!^)}  . 

(The  superscript  +  (or  -)  refers  to  the  saddle  point  at  cs  = 

(or  e,c  =  -Tf).) 

It  can  be  shown  that  for  <a~(g)  .>  10  Eq.  (30)  reduces  to 
Eq.  (27).  This  is  due  to  the  fact  that  F[Ka~(g)]  ->  1  as  ica±(g)  -v 
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but  for  practical  purposes  ica^e)  >  10  is  a  sufficient  condition 
to  replace  the  F  factor  by  unity. 

The  result  of  (28)  or  (29)  yields  a  uniformly  asymptotic  repre¬ 
sentation  for  the  total  field  so  that  the  diffracted  field  properly 
compensates  for  the  discontinuity  in  the  geometrical  optics  fields 
across  the  transition  boundaries,  thereby  yielding  a  total  field 
(geometrical  optics  plus  the  diffracted  fields)  which  is  continuous 
everywhere .  This  c?n  be  verified  by  taking  the  limit  of  Eq.  (28) 
or  (29)  as  one  approaches  the  transition  boundaries  from  either  side. 
Further,  the  total  field  obtained  via  the  limiting  process  on  either 
side  of  the  boundary  exactly  equals  the  "Cauchy  principal  value’1 
resulting  from  U(0)  =  ha  as  explained  earlier,  i.e.,  the  Cauchy 
principal  valui  yields  one  half  the  residue  contribution  associated 
with  the  pole  corresponding  to  geometrical  optics  incident  (nr 
reflected)  wave  when  it  happens  to  fall  exactly  on  the  saddle 
point  in  question  thereby  defining  the  total  field  on  the  shadow 
(or  reflection)  boundary.  This  continuity  of  fields  at  the  shadow 
and  reflection  boundaries  will  be  described  in  the  following  para¬ 
graph  . 

It  may  be  easily  verified  that  the  diffracted  field  just  within 
the  lit  side  (<j>-<j>'  =  u-e;  e>0)  of  the  shadow  boundary  gives  in  the 
limit  (e->0)  the  following: 


(34) 


/ 


v 


e-j[k(p+p')-ir/4]l 


» 


whereas  the  diffracted  field  just  within  the  dark  side  of  the  shadow 
boundary  {$-$'  =  ir+e,  e>0)  yields  in  the  limit  (e-+0)  the  following: 
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(35) 


y 


s-j  [k(p+p  ')-tt/4] 


The  incident  ray  field  which  is  visible  just  within  the  lit  side 
(<f>-<j>'  =  7t-g ;  as  e->0  and  e>0)  of  the  shadow  boundary  is 


(36) 


2 _  p-j[k(p+p')  -n/4  ] 

ihwt e 


(Equation  (36)  is  just  the  large  argument  approximation  for  the  line 
source  field,  -j/4  [k(p+p')] ) .  The  total  field  just  within  the 

lit  side  on  the  shadow  boundary  is  given  by  the  sum  of  the  results 
of  (34)  and  (36)  and  it  is  seen  to  equal  the  total  field  just  with¬ 
in  the  dark  side  of  the  shadow  boundary  given  by  (35).  Thus,  the 
total  field  is  continuous  across  the  shadow  boundary  and  equals  one 
half*the  geometrical  optics  incident  ray  field,  which  as  indicated 
earlier  is  exactly  the  Cauchy  principal  value  associated  with  the 
incident  wave  pole  for  the  shadow  boundary,  and  is  given  by 


i. 

f  i 

/  2  -j  rk(p+p  1  )-ir/4] 

2 

i 4 

J  jkwy  e 

* 

Similarly,  the  total  "scattered"  field  on  the  refaction  boundary 
is  one  half*the  geometrical  optics  field  of  the  r*,  icv-ced  ray  (just 
within  the  reflection  boundary).  The  incident  ray  field  is  of  course 
continuous  across  the  reflection  boundary.  The  total  field  (in¬ 
cident  plus  scattered)  is  thus  continuous  across  the  reflection 
boundary. 

It  is  of  interest  to  note  that  the  solution  presented  above 
for  the  canonical  problem  involving  a  perfectly  conducting  wedge 
excited  by  a  line  source  (electric  or  magnetic)  is  the  same  as 

*In  this  discussion  kpp'/p+p'  is  assumed  to  be  large. 
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ip 

that  conjectured  by  Rudduck  and  Tsai  except  that  the  results  of 
Reference  18  employ  the  Pauli  result  for  the  plane  wave  diffraction 
coefficient  instead  of  the  generalized  Pauli  result  given  in 
References  9  and  14. 


B.  Plane  Wave  Illumination 

The  plane  wave  incidence  case  is  treated  in  this  section.  The 
result  for  the  plane  wave  case  may  be  obtained  by  letting  p’  ->  »  in 
the  line  source  case  of  Section  A  so  that  kpp'/(p+p')  ■*  kp  in  the 
expression  for  k.  However,  we  shall  treat  this  problem  in  detail 
for  the  sake  of  completeness.  In  going  from  Eq.  (18a)  to  (18b)  the 
large  argument  approximation  for  KQ(jz(5))  was  used.  If  instead, 
an  approximation  corresponding  to  kp1  »  kp  is  used,  one  may  then 
wri  te 


K0(jk/p2+p,2-2pp'cos  zh  e  ^  ^  j^2irejkpC0S5  , 

as  given  in  Reference  4.  Thus  (18a)  becomes 

(38) 


[ji  ]£r  ^  ^  f 

If  one  suppresses  the  "line  source  factor  "  appearing  in  brackets, 
outside  the  integral  of  (38),  then  the  resulting  expression  which 
may  be  denoted  by  gp(p,B+)  is  used  to  obtain  total  field,  gp(p»<j>;<j>') 


due  to  a  plane  wave  normally  incident  on  the  edge  of  the  wedge. 

gj?  is  given  by 
h 


(39)  gP(p,W)  =  gp(p.6‘)  +  gp(p,6+), 

h 

where  for  a  plane  wave  of  unit  amplitude  and  zero  phase  at  the  edge, 

(40)  gW)i-V  f  artful  eik!>cos«d5. 

4njn  '  J 

L-L' 

Moreover,  gj?  satisfies  the  2-D,  homogeneous  scalar  Helmholtz  equation, 
h 

(41)  (v|  +  k2)  gP(p )  =  0, 

h 

the  Diriehlet  (s)  and  Neumann  (h)  boundary  conditions  on  the  surface 
of  the  wedge,  and  the  Meixner  edge  condition.  Thus,  for  an  electric 
plane  wave  field  zE  incident  on  the  wedge,  the  total  electric  field 
is  given  by 

(42a)  E  =  z  E1'  g$P  , 

A  ■? 

whereas  for  a  magnetic  plane  wave  field  zH  incident  on  the  wedge, 
the  total  magnetic  field  is 
(42b)  FT  =  z  H1  gjj. 

The  integral  appearing  in  (40)  has  been  extensively  treated  in 
Reference  14  where  the  generalized  Pauli  plane  wave  wedge  diffraction 

coefficient  was  obtained  via  the  method  outlined  in  Appendix  I.  One 
might  note  that  the  integral  over  L-L1  is  evaluated  exactly  as  done 
earlier  for  the  canonical  problem,  where  the  SDP (±tt)  closes  the 
L-L 1  contours  thereby  facilitating  the  use  of  the  Cauchy  Residue 
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Theorem.  Only  the  results  for  the  diffracted  fields  (SDP( ±tt)  con¬ 
tributions)  will  be  listed  here,  the  geometrical  optics  terms  being 

evaluated  via  the  Residues  of  the  integrand.  Thus,  the  diffracted 
od  .  . 

field  gj;  is  given  in  terms  of  the  scalar  plane  wave  diffraction 
h 

coefficient,  D$  as 
h 


(43) 


9S  '  D: 
h  I 


and  is  obtained  for  a  largness  parameter  k  =  kp. 


(44)  Ds  =  (d+($",n)  F[i<a+(ff)]  +  d“((f  ,n)  FCica“( B“) □)  + 
h 

+  {d+(s+,n)  F[Ka+(e+)]  +  d"U+,n)  F[ica"(6+)3}  , 


where 


-J 


d  ( B,n)  =  - 


e 


n/ZirrlT 


\  cot  (in^)  »  3  =  6+  5  <P  + 


Also, 


(45)  FOa^e)]  =  2jUKa±(e) 


J<a±(3) 


-4a±(s)j 


|  e”jT  dT, 


where  k  =  <p?and 

(46)  a±(e)  =  {1  +  cos(-p  +  2nN±iT)}  . 

The  transition  region  correction  factor  F[<a±(s)]  tends  to  unity  for 
<a±( 3)  >  10,  When  F  -*■  1  in  all  of  the  four  terms  in  (44),  reduces 
to  the  form  given  earlier  in  (27).  h 


One  obtains  (44)  by  utilizing  the  results  of  Appendix  I  for  evalu¬ 


ating 


L-V 


F1(c,8+)  eKf^ 


via  the  modified  steepest  descent  method,  where  F(s,B+)  and  f(^) 
are  taken  to  be 


1 _ 

4irjn 


f+B+  \ 

2n  /  J 


and 


f(s)  =  j  cos  ^,with  •<  =  kp. 


It  can  be  shown  that  the  complex  s-plane  topology  is  the  same  as  in 
Fig.  4  with  L-L '  and  SDP(±tt)  contours  unchanged.  The  integrand  of 
(40)  however  has  no  branch  point  singularities. 

C.  Spherical  Wave  Illumination 

Up  to  this  point,  only  2-D  geometries  have  been  treated.  Con¬ 
sider  now  the  case  of  a  scalar  point  source  illumination  of  the  wedge. 
Let  the  point  source  be  placed  at  s'(p' ,z').  The  point  source 
generates  scalar,  spherical  waves.  For  a  source  of  unit  strength, 
the  total  field  due  to  spherical  waves  incident  on  the  impenetrable 
wedge  is  denoted  by  gs(p,4>,z;p' ,4>* »z') ,  where  gs  satisfies 


(47)  (v2+k2)gs{i\?')  =  -  =-6(  js-'s"1 1 ) . 

h  p 


gs(s,s”')  is  the  scalar  point  source  Green's  function  for  the  wedge, 
h  _ 

relating  the  field  ats(p,4>,z)  due  to  a  unit  strength  point  source 
at  r(p*  ,z'). 
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gs(s,s  )  satisfies  the  Sommerfeld  radiation  condition,  the 
„  .  h 

Meix.ner  edge  condition  and  the  boundary  conditions: 

3gh(s,s'') 

*  0,  and  —  =  0. 

<£=0>nir  <f»=0  ,niT 


9s(s,s') 


Note  that  ?2  -  v2  +  t,  1s  the  3-D  Laplacian  operator.  This  3-D 
problem  may  be  reduced  to  a  2-D  problem  through  a  Fourier  integral 
transformation.  The  z -variation  may  be  removed  via  a  Fourier 
transform  on  the  z  variable  . 

Thus, 


(48) 


gs(p»$,h;p' ,<j>* ,z')  =1  g  dz 

h  i  c 


and  by  the  uniqueness  of  the  inverse  transformation , 

(49)  gs(s,s')  =  f  g$(p ><f>,h ;p *  ,<*>'  ,z')e^hz  dh. 
h  -«  h 

The  existence  of  the  transform  pair  is  assumed.  Thus,  Fourier 
transforming  (47)  yields 


(5°)  (V*  +  k£)J  (p,*,fi;p',**,2.) 

h 


=  --6-(p~p')4(<t-(l>f)  ~jhz' 


where 


k2  *  k2  -  h2. 


Rewriting  the  above  result  as 

(51)  {7?  +  kt%(p^;s,)ejhz'  = 

h 


-«(|F-P1), 
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% 


allows  one  to  interpret  gg  (p,<j>,h;s 1  )e^z  as  being  identical  to 

h 

( p ,p 1 )  (2-D  scalar  line  source  Green's  function  of  tq.  (4))  with 
tile  exception  that  k  must  now  be  replaced  by  kt  in  Eq.  (4).  Thus, 

gs  ejhz  may  be  replaced  by  Gs(^,p';  kt). 

h  h  _ 

The  integral  representation  for  G  (p,p‘ ,  k  )  is  known,  as  it 

is  equal  to  G  (p,p';k)  with  k  replaced  bv  k,  ,  and  is  given  by 

h  *  * 

(52)  Ss(P,^'  ;kt)  =  G(p,p‘  +  G(p,p'  *,8+;kt) 

h 

or 

V”2'  ’  9(6')eJhZ'  ?  <i(BVh2' 

h 

where 


g(6+)ejhz'  =  G(P,P';e+;kt)  =  [  |yCOt^jQ  K0(jz(0)<U  , 


as  given  by  (18a),  with  z (?)  being 


z(?)  =  kt/p2+p'2  -  2pp'  cos  K  =  X2-h2  Jphp'^-lpp 1  cos 


From  the  inverse  Fourier  transform  relation, 

(55)  g  (?,?')  =  g(p,z;p',z,;3")  +  g(p,z;p'  -’,*,(j'f) 
h  ’ 


where 


g(pvz;p' ,z';g) 


h  ff  |  TXT  «t(^)K0(oz(5,)e-^(2'-2>d51, 

;  U-L'  8n  Jn  J 
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(56) 


where 

Kv(jz) 


h(2) 


(z) 


has  been  used  to  obtain  (56)  from  the  preceeding  equation.  Further, 
19 

it  is  known  that 


dh  H^2)(/!?Ii?{p2+p,2-2pp'  cos  ah)  e~ih{z'-z) 
-ikv42+p,2-2pp'cos  5+(z-z‘)2 


=  2j 


A p 2  +  p‘2  -  2pp‘  cos  5  +  (z-z1)2 


Ut’'1'izing  the  above  result  in  (56)  yields: 


(57) 

§( P  s  Z  jp  ,Z  5S) 


e- J  k>4  2+p 1 2-2pp 1 cos£+ ( z-z 1 ) 2 
v^p+p  1 2-2pp 1  cosE+( z-z 1 ) 2 


The  usefulness  of  the  above  point  source  scalar  problem  in 
the  treatment  of  the  vector  electromagnetic  problem  is  seen  if 
one  considers  the  case  of  either  an  electric  or  magnetic  dipole 
moment  (e+Jut  time  dependence  assumed)  which  is  z  <  rected,  and 
radiates  in  the  presence  of  a  wedge.  The  magnetic  vector  potential, 
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ft  at  s  (observation  point)  for  an  electric  dipole  moment  of  strength, 

A 

zpg6( i s  -s'|)  located  at  s',  is  given  by 

(58)  7(7)  =  zupe  gs(s,‘s') 

whereas  the  electric  vector  potential,  ?  at  7  for  a  magnetic  dipole 

A 

moment  of  strength,  zpm<S(|s  -s'|),  is  given  by 

(59)  ?'(?)  =  zepm  gh(?,I'). 


Also, 

(60) 


TT 


VXA 

u 


and  E„ 
m 


> 


where  the  subscripts  e  and  m  are  used  to  denote  field  quan  ties 
resulting  from  an  electric  dipole  moment  or  a  magnetic  dipole 
moment,  respectively.  Thus,  =  magnetic  field  intensity  due  to 

A 

zpg6(|s  -s' | ) ,  whereas  ^  =  electric  field  intensity  due  to 

A 

zpm6 ( j s  -s' |).  An  arbitrarily-polarized  spherical  wave  can  be 
generated  by  a  superposition  of  the  fields  of  these  two  dipoles. 

Figure  5  depicts  the  geometry  involving  a  point  source  at 
7'  and  the  observation  point  at  7,  for  the  source  radiating  in  the 
presence  of  the  wedge  which  is  impenetrable.  A  coordinate  transforma¬ 
tion  is  now  introduced.  One  is  referred  to  Fig.  5  for  details. 

Let  p'  =  s'  sin  bq  and  p  =  s  sin  Bn.  Hence,  j(z-z')j=(s’+s  )cos  BQ, 
and 


(61) 


)2+D'2-2pP'cos  5  +  (z-z1)2  =  (s'+s  )2  sin2  bq  - 

-2s' s  sin2e0(l+cos  ?)+(s'+s  )2cos2 


=  (s'+s  )2 


2s' s  sin2B 


r  zs  s  s 

L1-  — 


- -  (1  +  COS  c)|  . 

s  )2 


i 
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Incorporating  the  above  result  in  (57)  yields 


(62)  g(p,z;p',z';B)  =  — ^ —  [  dg  cot  ■ 

16-rr^jn  JL-L'  \  ZnJ 


- * - 

2s' s  sin^B 

-jk(s*  +s  )  /  1 - - - (1  +  cos  t) 

U  +s  r 


J  (s'  +s  )2-2  s'  s  sin2e  (1  +  cos  z) 


As  before,  the  integral  over  L-  L'  is  evaluated  asymptotically  for 
s'  s  2 

large  k  sin^,  and  the  SDP(±tt)  again  close  the  L-L'  contours 

allowing  the  use  of  the  Cauchy  Residue  Theorem.  The  location  of  the 
SDP  through  saddle  points  at  Zc  =  ±ir  and  the  complex  z  plane  topology 
is  once  again  unchanged  except  that  the  branch  points  of  the  integrand 
are  now  given  by 

,  ^2+s2  +  2s‘  s  cosV 

Zh  =  2u  ±  j  cosh"  - - - 2 - -  >  *  =  0,±1,±2»* •  • 

2s  s  sin  3 


and  one  is  referred  to  Fig.  4  for  details. 

The  asymptotic  evaluation  of  the  integral  as  in  (62)  is  per¬ 
formed  via  the  modified  steepest  descent  method  outlined  in  Appendix 
I,  withF^)  and  f(^)  now  given  by 


F,U)  =  -^V~ 
1  16-n^jn 


cot(  S6) 


-jk(s‘  +  s  ) 

e _ _ 

/(s'  +s  )-2s' s  sin2B0(l+cos  5) 


9 


and 


f(z)  =  J  [1  +  cos  z] 
where  the  approximation 
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(1  +  cos  c) 


* 


-jk(s'  +s  ) 


1  - 


(63) 


2s' s  sin26c 
(s  +s  )2 


-jk(s  +s  )  ^  1  - 


«  e 


i  •  2n 

s  s  sin  6, 


(s'  +  s  )2 


(1  +  cos  e)r 


has  been  used  in  (62)  for  reasons  identical  to  those  involved  in  the 
approximation  of  Eq,  (19). 


(  SIDE  VIEW  ) 


Qe  IS  A  UNIQUE  POINT  ON  THE 
EDGE  FOR  A  GIVEN  SOURCE 
AND  OBSERVATION  POINT. 


Fig.  4,  Geometry  for  a  point  source  illumination 
of  the  wedge. 
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The  residues  corresponding  to  the  poles  of  F^U)  at  *■  =  £p  = 
-0  +  2nNn  (N-0,  ±1,  ±2,“*)»  yield  the  geometrical  optics  fields 
(Note  that  |cp|  <_u).  The  saddle  point  results  yield  the  dif¬ 
fracted  ray  field  component. 

As  before,  the  diffracted  fields  may  be  denoted  by  gd(?,?'), 

h 

where 


(64)  gd(s,s')  =  gd(p,z;p'  .z'je")  +  gd(p,z*,p'  ,z'  ;$+), 
h 


with 


gd(p,z;p' ,z' ;e?)  =  -  j  FjU,®*)^ ^  -  F1(§,0+)eKf . 


SDP(tt) 


SDP(*) 


s  s 


7 

(Note  that  k  =  k  ^r^r-  sin  in  tiis  case). 


No  derivations  are  presented,  but  it  can  be  shown  by  using  the  results 


of  Appendix  I,  that 


-jk(s  +s  ) 


SDP(±tt) 


2tt(  S*  +S  ) 


(  TT 

V 

nr 


16ir4jn(s  +s  )  J  ks  s  sin2e. 


•  cot  (^)  pC<a±lB)D 


which  may  be  re-written  as: 
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(65)  -  J  F1(5,e)eKf(e)d^  ij^r 
SDP(±u)  L 


:[Ka-(s)] 


,-j  f/4 


I  cot  (4  iT? 


One  is  now  able  to  interpret  the  asymptotic  result  in  terms  of  the 
GTD,  through  a  diffraction  coefficient,  Ds.  Thus 


(66)  gf  'v  u1  -  ~ 


s  s'  (s'  +s  ) 


v/here 


u1  =  field  of  the  unit  source  incident  at  the  edge,  and  is 


given  by 


.1  e 


,  (and  is  equivalent  to  the  point  source,  scalar 
Green's  function  in  the  absence  of  the  wedge). 


(67)  Ds  =  [{d  (6",n)FC«a  (ff)3  +  d“(g  ,n)F[Ka  (g  )]}  + 
h 

?{d+(6+,n)F[Ka+(8+)]  +  d"(e+,n)F[:Ka"(6+)3}]  , 

sin  8q 


(68)  d  fg,n)= 


icot 

n/S  2  2"  ’ 


B=e  =  <t>  + 


FCtca'^e)]  =  ejKa"(3)  f  e‘jT  dx 

i^^gj 
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in  which 

(69)  a±(e)  =  {1  +  cos (-6  +  2nN±Tr)}>  k  =  sin2  3Q, 

and  N±  denotes  the  poles  closest  to  each  of  the  two  saddle  points. 
Once  again,  F  may  be  looked  upon  as  a  transition  region  correction 
factor  which  tends  to  unity  away  from  the  transition  boundaries  or 
whenever  k  a±(e)  >  10.  The  geometrical  optics  contribution  may 
easily  be  shown  to  comprise  of  terms  of  the  type 


-jk/p2+p,2-2pplcos[-3+2nN-TO  + 

(70)  1. - _ , - - —  — - -  U[ir-|-6+2nN±Tr|  ]. 


4 


4tt|o2+p  1 2-2pp 1  cos[-3+2nN“nI] 


If  F  *>  1  in  (67),  then  Dg  of  (67)  reduces  to  Ds  as  given  in  (27). 


The  (sin  B0)_1  factor  associated  with  Ds  in  (67)  indicates  a 

h 

conical  spreading  of  the  rays  in  space  as  schematized  in  Fig.  5. 

This  conical  spreading  of  the  diffracted  rays  is  a  consequence  of 
Fermat's  principle. 

Summarizing  the  results  derived  thus  far,  one  observes  that 
the  canonical  problem  of  a  line  source  wedge  illumination 
(cylindrical  wav.  incidence)  is  useful  in  treating  the  subsequent 
scalar  problems  involving  plane  wave  and  spherical  wave  illuminations 
of  the  wedge.  Separating  the  line  source  factor  out  from  the  in¬ 
tegral  formulation  of  the  canonical  problem  yieldsthe  solution 
for  the  plane  wave  incidence  case,  whereas  a  Fourier  transformation 
related  the  spherical  wave  incidence  case  to  the  cylindrical  wave 
incidence  case. 
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The  asymptotic  approximation  in  each  of  the  problems  treated 
results  in  a  solution  which  can  be  interpreted  in  terms  of  the 
geometrical  theory  of  diffraction,  and  as  expected,  outside  of  the 
transition  regions  of  the  shadow  and  reflection  boundaries,  the 
diffraction  coefficient  is  independent  of  the  type  of  edge  il¬ 
lumination.  Within  the  transition  regions,  where  strictly  speaking, 
ray  optics  is  no  longer  valid,  the  geometrical  theory  of  diffraction 
can  be  used  formally  if  one  corrects  the  diffraction  coefficients  in 
the  following  simple  manner.  The  scalar  diffraction  coefficients 
may  be  written  as  the  sum  of  four  terms,  and  each  of  these  terms 
is  multiplied  by  a  function  F[xa],  see  Eq.  (32)  for  example. 

|F(Kd)|  <  1; 

outside  of  the  transition  regions,  where  xa  >  10.  F[xaH  ^  1. 

It  is  very  interesting  that  FC<aD  has  precisely  the  same  form 
for  the  different  edge  illuminations.  The  only  difference  for  the 
different  types  of  illumination  occurs  in  it's  argument;  if  we 
write  xa  =  kLa,  then 

p  for  plane  wave  illumination 


L  = 


pp 


-r-  for  cylindrical  wave  illumination 


s  +  sr  sin  8o  ^0r  sP^encal  wavs  illumination. 

To  ensure  the  validity  of  the  asymptotic  solution,  kL  should  be  large. 

It  is  indeed  remarkable  that  the  complex  field  behavior  in  the 
transition  regions  can  be  described  in  this  compact,  simple  way. 
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In  the  next  chapter  we  shall  show  how  these  scalar  (acoustic) 
problems  are  related  to  the  electromagnetic  problem  of  a  perfectly- 
conducting  wedge  illuminated  by  an  arbitrarily  polarized,  plane, 
conical,  cylindrical  or  spherical  wave.  It  will  be  seen  that  the 
scalar  diffraction  coefficients  developed  in  this  chapter  play  an 
important  role  in  the  dyadic  diffraction  coefficient  of  the 
electromagnetic  problem. 
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CHAPTER  III 

ELECTROMAGNETIC  WEDGE  DIFFRACTION 

In  Chapter  II,  the  plane  waves  and  the  cylindrical  waves 
associated  with  the  uniform  line  source  are  normally  incident 
on  the  edge  of  the  wedge.  The  solutions  of  these  scalar  problems 
can  be  directly  related  to  the  corresponding  electromagnetic  prob¬ 
lems,  see  (5)  and  (6)  for  the  line  source  case  and  (42a, b)  for  the 
plane  wave  case.  In  this  chapter  the  normal ly-inci dent  waves  of 
Chapter  II  will  be  generalised  to  obliquely-incident  electro¬ 
magnetic  plane  waves  and  obliquely-incident  conical  electromag¬ 
netic  v/aves  whose  sources  are  traveling-wave  line  currents. 

In  addition,  the  scalar  spherical  wave  excitation  will  be  gen¬ 
eralized  to  an  electromagnetic  spherical  wave.  Arbitrarily- 
polarized,  plane,  conical  and  spherical  waves  may  be  treated  by 
properly  superimposing  solutions  of  the  electric  and  magnetic 
type.  Asymptotic  solutions  to  these  more  general  problems  and 
their  interpretation  in  terms  of  the  geometrical  theory  of  dif¬ 
fraction  form  the  material  of  this  chapter. 

However  before  taking  up  these  problems,  a  ray-fixed  co¬ 
ordinate  system  will  be  introduced.  It  will  be  seen  that  when 
this  coordinate  system  is  used,  the  resulting  dyadic  diffraction 
coefficient  may  be  written  compactly  as  the  sum  of  two  dyads, 
or  alternatively,  it  may  be  written  in  matrix  notation  as  a 
2x2  diagonal  matrix. 
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A,  Coordinate  Systems 

The  circular  cylindrical  coordinate  system  (p,<bz;p'  *<J>'  ,z‘) 
which  we  have  useu  to  treat  wedge  diffraction  is  an  edge-fixed 
coordinate  system.  On  the  other  hand,  the  rays  incident  and 
diffracted  at  in  Fig.  4  are  more  conveniently  described  in  terms 
of  spherical  coordinates  centered  at  Q^.  Let  the  position  of  the 
source  of  the  incident  ray  be  defined  by  the  coordinates  (s',g^,<j>') 
and  the  point  of  diffraction,  by  the  coordinates  (s,B0,<j>)  as  shown 
in  Fig.  5.  The  orthogonal  unit  vectors  associated  with  these 

A  A  /.  AAA 

coordinates  (s *  #3^ 1 »  s,B0,<j>)  are  ray  fixed  in  contrast  with  the 

A  A  A  AAA 

unit  vectors  (p1  t<t>‘  ,z'  ;p ,<j>,z). 

The  plane  containing  the  incident  ray  and  the  edge  of  the 
wedge  will  be  referred  to  as  the  plane  of  incidence,  and  the  plane 
containing  the  diffracted  ray  and  the  edge  of  the  wedge,  as  the 

A  A 

plane  of  diffraction.  The  unit  vector  I  =  -s’  is  in  the  direction 
of  incidence  and  the  unit  vector  s  is  in  the  direction  of  di f- 

A  A 

fraction.  It  is  apparent  that  the  unit  vectors  6^  and  are 
parallel  and  perpendicular,  respectively,  to  the  plane  of  inci- 

A  A 

dence,  and  the  unit  vectors  and  $  are  parallel  and  perpendicular, 

A  A  A 

respectively,  to  the  plane  of  diffraction;  moreover,  =  <f>'  x  I 

A  A  A 

and  gQ  =  <j>  x  s.  It  will  be  shown  later  that  when  the  incident 
field  is  resolved  into  its  parallel  and  perpendicular  components 
and  the  diffracted  field  is  resolved  into  its  parallel  and  per¬ 
pendicular  components,  the  diffraction  coefficient  can  be  ex¬ 
pressed  as  a  2  x  2  diagonal  matrix,  and  in  this  sense  the  ray- 
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fixed  coordinate  system  is  a  preferred  coordinate  system. 

In  the  following  sections  the  diffracted  field  components 
will  first  be  determined  in  the  edge-fixed  coordinate  system 
and  then  transformed  to  the  ray-fixed  coordinate  system. 

B.  Plane  Wave  Illumination 

The  case  of  an  arbitrarily-polarized,  electromagnetic  plane 
wave  obliquely  incident  on  a  perfectly-conducting  wedge  is  treated 
first.  The  geometry  of  the  problem  is  depicted  in  Fig.  5.  The 
incident  electric  field  F1  may  be  decomposed  into  a  component 
parallel  to  the  plane  of  incidence  E„  =&^  *  F  ,  and  a  component 

perpendicular  to  the  plane  of  incidence  Ex  =  4. '  •  F  .  Thus, 

F7  =  0^  E,^  +  4> '  Ejj,  and  H7  =  0^  H7  +  4. 1  The  z  componentsof 
E,7  and  EjJ  are 

(71)  E(]z  =  z  •  F7  =  ($£  .  F7)  sin  0Q  and  E^z  =  0. 


A 

The  z  components  of  the  magnetic  fields  associated  with  the 
parallel  and  perpendicular  components  of  the  incident  electric 
field  are 

(72)  H7,  =  0  and  H1’  =  M  (J*  •  E7)  sin  a  . 

11 Z  1 Z  1  e  0 

(Note:  B7  =  J'  •  F7  and  H]_  =  3^  •  K7) 

The  total  electric  field  intensity  F  and  the  total  magnetic 
field  intensity  IT  in  the  region  surrounding  the  wedge  may  be  de- 

A  ^ 

composed  into  total  transverse  (to  z)  fields  O^JT^.)  and  total 
axial  fields  (Ez>Hz),  respectively.  Thus, 
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* 


(73)  IT  -  +  z  Ez  and  IT  =  FTt  +  z  Hz. 

Let  GS(p»4>,z;<f>')  represent  the  total  field  due  to  a  "scalar" 
h 

plane  wave  obliquely  incident  on  an  impenetrable  wedge,  G°(p,<j> jZj^' ) 

n 

satisfies 

(v2  +  k2)  G°(p,*fz;4>')  =  0. 
h 


A  traveling  wave  dependence  in  z  allows  one  to  express  G°(p;<j>')  as 

f  -jk  z  5 

G°(p;<}>')  =  gj?(p;4>’;kt)  e  1  ,  where 

h  n 


(74)yz  51  k  cos  eo,  (v2  +  ky  g^(p ,4>;4>;kt) 

h 

and  k^  =  k  sin  6Q. 


0, 


9P(p ;4> '  . )  is  identical  to  the  expression  for 

h 

in  (39)  with  k  replaced  by  kt. 

It  follows  that 


f 


(75)  ^and 


■  Hz  4 


If  one  knows  the  solutions  to  the  axial  fields  Ez  and  Hz, 
then  one  can  evaluate  and  H"t  in  terms  of  the  axial  fields  from 
the  results  of  Appendix  II,  given  by  the  folowing  relations 


v.  E  z  x  V.H 

ci,  t  Z  ±  c  uZ 

•Jkz  7TJ  +  3w 


(76) 


k  -k; 


k2-k2  ’ 


V 


Vt  H  Z  X  VtE 

-jk,  4-4  -3“s 


'2  k2-k2 


k2  -  k2 
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The  above  expressions  are  valid  for  problems  which  possess  cy¬ 
lindrical  uniformity  in  structure  and  fields. 


The  solutions  for  Ez  and  Hz  are  already  known  from  Chapter  II. 

The  z  components  of  the  diffracted  fields  Ed  and  Hd,  may  be 

written  in  terms  of  the  diffraction  coefficient  D  of  (44),  but 

n 

with  k  replaced  by  k^.  The  total  axial  field  is  of  course  a 
superposition  of  the  geometrical  optics  field  and  the  diffracted  ray 
field.  Since  this  report  is  primarily  concerned  with  the  calculation 
of  diffracted  fields,  only  the  results  for  the  diffracted  field  will 
be  presented  here.  It  is  a  straight-forward  matter  to  obtain  the 
geometrical  optics  contributions  via  the  residues  at  the  poles  as¬ 
sociated  with  the  incident  and  reflected  fields,  as  was  done  in 
Chapter  II,  or,  alternatively,  the  geometrical  optics  contributions 
can  be  deduced  directly  from  physical  considerations. 

For  the  diffracted  ray  contribution  to  the  axial  fields, 
gj?(p»<f>;4>'  ;kt)  is  replaced  by  gPd  as  in  (43)  (but  with  k  replaced 

by  k.).  gPd  is  in  turn  expressed  in  terms  of  D  of  (44),  but  with 

n  h 

k  replaced  by  k,.  ihus, 

-d  .  ri  fTTTT  n  e  r  .  JKz2 

9 


(77) 


Ez  *  E\  Jsir\$ 0  D$  -j=-  e 

and  -jk+p 

Hd  *  HMsin  0“  D 


tk  -ok  z 
e 


z  z  v  -  ^o  h  / — 

>1  p 

where  E1 ,  H1  are  evaluated  at  Qr,  the  point  of  incidence  on  the  edge, 

Z  Z  £ 

(78)  Ds  =  sln~V~  C{d+(3"»n)F[;Ka+(e)]  +  d” ( ,n) FCo" ( 3")  + 

h  °  +  (d+(g+,n)F[<a+(B+)3  +  d“(g+,n)F[Ka"(8+)3)3  , 
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* 


with  k  =  ktp 
d±(s,n) 


kp  sin  e  , 

_e"J  f 


n fib 


1  cot  ( !=§-' 

2  C0T:  V  2n 


F[Ka*(B)]  =  2j|Jicai(B)'|  ej*a±(f3) 


fjT  dv 


lJKa£(sT| 


and  a“(g)  =  1  +  cos{-8  +  2nN~n),  8  =  3“  =  <(>  ±  <j>' .  One  may  refer  to 
Chapter  II  for  the  definition  of  N*.  From  Fig.  5  it  is  evident  that 
p  =  s  sin  8  .  Also, 

e-j(ktp+kzz)  _  e-jks  (w-j th  as  the  origin). 


_  _  * 

To  evaluate  E^.  and  one  may  take  the  vector  and  z  x  v.  operators 

within  the  integral 

r  3\p  cos  5 


9P(p,*iV;kt)  =  5^. 


cot 


L-V 


(H^J 


d? 


where  gp(p ;kt)  =  gp(p,$;*'  ;kt)  +  gp(P, ;kt) 
h 


(Note:  as  before,  gp 


=  0,  and 


Thus, 


(79)  { 


cj!=0,mi 

ci  -sv 

V  E  %  o  -—r — , - 1  (  i 

t  z  4ifjn 


£ 
3  4> 


J 


■  ^  j  jkt  cos  f,  e 

•[cot(Hr-)-  cot(Hr-)]dt 


=  0). 

<{>=0,mr 

jktp  cos  e 


i  -jkzz  ( 
Hi  e  z  [ 


_!!*»»  A  2  i  (  St. 

Vt"z%0~~j J  JKtwa* 

[L-f 

•H^)  +  cot(4^)]deJ, 


jktp  cos  ^ 
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»  •  •  • 
where  Ez  and  are,  respectively,  E,1,  sinpo  and  sinj30  taken  at  Q^. 

The  approximation  involved  in  (79)  is  such  that  only  those  terms 
are  retained,  which  on  asymptotic  evaluation  yield  terns  of  0^  . 

Higher  order  range  dependent  terms  are  neglected;, this  is  a  valid 
approximation  if  kp  is  sufficiently  large.  The  diffracted  field  con¬ 
tribution  to  v^Ez  and  vtHz  (denoted  by  v^E^  and  v^Hz  respectively)  is 
obtained  from  the  asymptotic  approximation  of  (79)  via  the  modified 
method  of  steepest  descent  described  in  Appendix  I;  G.oe  finds  that 


(80) 


vj: 


■d 


Vzj 


H 


tzU-Ph-jk  sin3/2  .)  D  a 


-jkz 


F 


v/ith  D  in  the  above  equation  being  identical  to  that  of  D  in  (78). 

«  h 

One  may  re-write  (77)  and  (80)  in  terms  of  the  parallel  and 
perpendicular  components  of  the  incident  field  as 


(81) 


r  ,  ^ 

r  F 

u 

E1 

2 

n  p 

”z 

\  4  J 

X, 

*  J7 


(82) 


(  a' 

f  .  \ 

E1 

ii 

r 

7tEz 

1 

* 

L~p  J 

,’tHz, 

Hi 

iHi  J 

iJl, 


U-jk  sin  30)  D, 


-jks 


F 


Utilizing  the  above  equations  in  (76),  noting  that  p  =  s  sin  $0  + 

A  A  A  A  AAA 

3Q  cos  S0,  7.  =  s  cos  -  e0  sin  Bfl  and  z  x  p  =  <>,  one  obtains  the 
following  expressions  for  the  diffracted  fields  (E^Ji^)  expressed 
in  terms  of  the  coordinate  system  fixed  in  the  diffracted  ray 

(33)  E*  .  +  itl  •«.  -60  E,;  Ds 


e"jks  *  • 

"  ♦  Ex  Dh 


•jks 


F 


F 
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(84)  «-  -  R?  -W  Oh  ^  -  5.  H’  Os  ^ 

>/  s  JS 

The  relations  h|  J"u/e  =  -  F.|  and  hJ  Jv/e  =  E.,  have  been  used 
to  express  the  above  results  for  and  in  terms  of  (eJ,E^)  and 
)  respectively. 

In  equations  (83)  and  (84)  the  incident-field  is  resolved  into 
components  parallel  and  perpendicular  to  the  plane  of  incidence.  Now 
if  the  diffracted  field  is  expressed  in  terms  of  its  components  parallel 
and  perpendicular  to  the  plane  of  diffraction, 

t3  =  e$0  +  e?;  , 

^  ■  H$o  +  H«*  - 

one  may  rewrite  the  above  expressions  in  the  form  of  the  geometrical 
theory  of  diffraction. 


and 


.  -  _  ~jks 

E'fQr)  •  Dp  - - 

E  E  JJ- 
-jks 


•»  hVQj)  •  dh  £ 


JT 


where  the  dyadic  diffraction  coefficients  Dc(  4>»  <p' ;  8  )  and  Du(  <f>,  <J>' ;  3  ) 

tOHO 

are  expressed  in  the  following  simple  form 


De(*.*'.B0)  =”®o®o  D  ( 4> , «f> 1 ,3  )  -<J> * Dl. ( 1 >8r ) 

(85)  =  0  °°  S  _ 

>80)  ="^o®o  *tf0) 
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where  Ds(<m',bo)  is  given  by  (78). 

Expressing  the  diffraction  coefficients  in  matrix  notation 


0  ' 

rEj(QE>' 

e”jks 

. 0  -V 

f 

K 

\  0  ' 

’Hi<  qe>’ 

e-jks 

k 

,°  -°s. 

,HX>. 

The  diffraction  coefficient  reduces  to  a  2x2  diagonal 
matrix,  because  we  have  expressed  the  field  components  in  the  proper 
coordinate  system,  namely,  the  ray-fixed  coordinate  system,  and 
because  we  retained  only  the  range  dependent  terns  of  0(  1/,/p). 

Keller  and  Lewis  have  obtained  the  dyadic  diffraction  coefficient  for 
a  wedge  in  an  edge-fixed  coordinate  system,  where  the  incident  and 
diffracted  fields  are  expressed  in  terms  of  their  cartesian  com¬ 
ponents.  They  obtain  a  3  x  3  matrix  with  7  non-vanishing  elements. 
Furthermore,  unlike  the  results  given  here,  their  diffraction  coeffi¬ 
cient  is  not  valid  in  the  transition  regions  at  the  shadow  and  re¬ 
flection  boundaries.  We  have  transformed  our  dyadic  diffraction 
coefficient  to  the  edge-fixed  cartesian  coordinate  system.  Outside 
of  the  transition  regions,  where  the  correction  factors  may  be 
replaced  by  unity,  our  dyadic  diffraction  coefficient  is  identical 
with  that  of  Keller  and  lewis,  which  shows  that  the  degree  of  ap¬ 
proximation  in  the  two  asymptotic  solutions  is  the  same.  In¬ 
dependently  of  our  work,  L.J.  Kaplan*  has  shown  that  the  Keller 

*L.  J,  Kaplan,  private  communication. 
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and  Lewis  dif fraction  matrix  can  be  reduced  and  diagonalized  by  a 
suitable  matrix  transformation. 

C.  Conical  Wave  Illumination 

Consider  next  the  diffraction  of  a  conical  electromagnetic 
wave  by  a  perfectly-conducting  wedge.  At  the  edge  the  direction 
of  incidence,  which  is  normal  to  the  conical  wavefronts,  is 
given  by  the  unit  vector  I;  I  •  z  =  cos  8  .  where  -rr/2  -  is  the 
half-angle  of  the  conical  wavefronts.  Conical  waves  may  be  gen¬ 
erated  by  the  electric  and  magnetic  traveling-wave  line  currents 
a  -jk  z'  a  -jk  z' 

zl  e  *  and  zM  e  ,  respectively,  where  I  and  M  are  complex 
constants  and  k  =kcos  8„. 

-Jk  z 

The  axial  (z-directed)  incident  electric  field  due  to  le 

-jk  z 

and  the  axial  incident  magnetic  field  due  to  He  may  be  shown 
to  be. 


(87) 


-jk,z' 
Ie  z 


-jk  z’ 
Me  z 


[f  Ho2)(V)] 


,  and 


,  respectively 


with 


k^  =  k2  -  kl  =  k  sin  eQ. 


The  primes  are  used  to  denote  a  source-fixed  coordinate  system. 

Next  let  us  examine  the  conical  wave  behavior;  the  electric 
current  line  source  will  be  treated.  The  fields  of  this  source  are 

O 

k  *  p  ,  . 

1  —  /  o  \ 

H-I(ktp,)| 


..  k! 

F 1  ss 

Z  Jtoe 


(88) 


le  ‘ 


=  0 


pi  _  *i  ^  ^  ^i  _  %  jb)£  <>  e 

Et  "  p  .2  3p'  *  9  t  *  t2  3^rtz* 

Kt  t 
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For  ktp'  large,  the  Hankel  function  may  be  replaced  by  its  asymptotic 
approximation,  so  that 


ktp*  +  ~  ^p'  eo  +  z*  coc;  6o^  =  constant  describes  a  conical 

surface  and  I  •  E"  and  I  •  FT  both  vanish;  consequently,  for  k^p ’ 
large  enough  for  the  Hankel  function  of  order  zero  to  be  replaced 
by  it's  asymptotic  approximation,  the  electromagnetic  field  of  a 
traveling-wave  electric  current  line  source  is  a  conical  wave  with 
it's  field  vectors  tangent  to  the  wavefront,  i.e„,  perpendicular 
to  the  associated  rays.  The  electromagnetic  field  of  the  traveling- 
wave  magnetic  current  line  source  has  the  same  properties.  In  the 
case  of  the  electric  current  line  source,  the  incident  electric 
field  lies  in  the  plane  of  incidence,  whereas  in  the  case  of  the 
magnetic  current  line  source,  it  is  perpendicular  to  the  plane  of 
incidence. 

The  total  axial  fields  (Ez  and  Hz)  excited  by  electric  and 
magnetic  traveling-wave  line  sources  can  be  shown  (see  Appendix 
III)  to  be  0 

r  -jk  Z 

(90)  Ez  =  Ie  2  and 

< 

k2  —jk  2 

H  =  -4-  Me  z  G.  (p",p*'  ;k. ) ,  respectively. 

^  Z  JWU  H  t 
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The  special  Green's  function  Gs(p,p'*,kt)  is  the  same  as  that  of  (4) 
with  k  replaced  by  kt#  As  before, 

(91)  Gg(p,p'';kt)  =  G(p,p' ;g";kt)  +  G(p,3';  +;kt), 
h 

ft 

where  the  integral  representation  of  6(p,p';B  ;kt)  is  given  by  (18b) 
with  kt  used  in  place  of  k.  Also,  the  approximation  of  (13)  is 
utilized  in  the  exponent  of  the  integrand  in  (18b).  Specifically, 


jkt(p2+p' 

(1  +  cos 

d£  . 

The  diffracted  ray  contribution  to  (91)  is  known  from  the  result 
of  (28)  in  Chapter  II  from  which  the  axial  diffracted  fields  may  be 
written  as 


where  D_  is  identical  to  that  of  (78),  but  with  k  equal  to 


(92)  G(p,p' ;B+;k.)  ^  ^ —  [f 

8tt  jn  /.-{.n  I 

-jk.(p+p'){l  -  ^ 

.  „  t  (p+p *)2 


k.  pp ' 


KPP 
P+P  ' 


sin  3 


o* 


In  terms  of  E1*  and  which  are  the  fields  incident  on  the  edge  at  QF, 
z  z 

it  is  easily  seen  that  (93)  can  be  written  as 
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(94) 


•jktP 


■j'M 


s 

h 


F 


JTuTc 


The  transverse  components  of  the  diffracted  fields  may  be  ou¬ 
tlined  via  (76)  in  the  manner  similar  to  that  of  the  plane  wave  case 
treated  earlier.  The  vector  and  z  x  operators  are  taken  within 
the  integral  of  (92)  before  carrying  out  the  asymptotic  evaluation. 

K4 '  4 »,  *^0  ^  ■ 


(95H 


l'tHz  -  -i*  4  Di>  s™3'\ 1 


"'*ktp  ~jk7z 
e 


F 


where  Dg  is  identical  to  that  of  (78)  but  with  <  -  sin  bq. 

h  _ 

In  deriving  (95),  terms  of  0(l/J~p)  are  retained;  all  higher 

order  range  dependent  terms  are  negligible  for  kfcp  sufficiently  large. 

One  might  observe  that  (94)  and  (77'  possess  identical  forms  as  do 

(95)  and  (80),  respectively. 

Since  =  E,?  and  in  this  problem,  =  E^  sin  gQ  = 

E1'  sin  6  .  and  =  H?  sift  e_.  Also,  e]  =  -/y/e  ,  and 

1  °  ,  z  "  °  -j(k.p+k  z) 

Ejj  =/ y/e  !■{, .  Utilizing  these  results  and  the  fact  that  e 

e~jks,  (see  Fig.  5),  one  may  express  (94)  and  (95)  respectively  as 


(96) 


k] 

< 

■< 

*  1 

Hd 

'V/ 

H/ 

l  ZJ 

\  ~ 

Ds  sin  6o 


-jks 


js 


,  and 


(Note:  the  subscripts  e  and  m  on  E^  and  respectively  denote  the  fields 
of  electric  and  magnetic  traveling -wave  line  sources.) 
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r  w*1 

\<) 

(97)  < 

v 

^  P  (“jk) 

H.1 

^  t  2, 

s 

h 


sin^  B 


,-jks 


0  jr 


Equations  (96)  and  (97)  are  similar  in  form  to  (81)  and  (82), 
respectively. 

Following  the  procedure  indicated  for  the  vector  plane  wave  in- 
cidence  case,  one  may  utilize  (96)  and  (97)  in  (76)  to  obtain  the 
total  diffracted  fields  in  the  ray-fixed  coordinate  system.  Then 
(I?**  and  fl^)  may  be  written  in  the  form  of  the  geometrical  theory  of 
diffraction  as 


_  -jks 

Z*  ^e\Qe)  *  De  £pr-  ,  and 
J  s 

*  F1  (Qp)  •  F  6  J 


H 


JT  * 


Here,  0^  are  the  dyadic  wedge  diffraction  coefficients  for  a  conical 
K 

wave  excitation. 


\  *  Ds  -  °h  > 
H  h  s 


and 

D  =  •rr^-HI(d+(g",ri)rLKa+(B'<)]  +  d"(B”,n)F[;Ka  (&  )3)  + 

l  b1n  8o 

+  {d+(8+,n)FCKa+(sT)3  +  d  (B"r,n(FCiea  (bt)3)D  , 

where  the  d±(B,n)  terms  are  identical  to  the  d±(B,n)  given  in  (78) 
for  the  plane  wave  case.  The  form  of  F{jca“(B)3  is  the  same  as  in 
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(78)  for  the  plane  wave  incidence  case,  except  that  k  in  the  cor¬ 
rection  factor  F[>al]  is  equal  to  sin  ba  in  this  case.  As  the 

p-r-p  o 

line  source  is  removed  from  the  edge  p'  -*•  »,  the  local  illumination 
is  plane  wave  and  <  -►  k  sin  bq,  as  it  should. 

When  eo  =  tt/2.  which  corresponds  to  uniform  line  source  ex¬ 
citation  (cylindrical  wave  incidence),  then  D$  are  given  by  (30). 
Finally,  as  the  arguments  of  the  transition  region  correction  factors 
become  large,  i.e.,  when  <a  >  10,  the  diffraction  coefficients  are  equal 
to  the  given  by  (27)  multiplied  by  (sin  80)~*. 

The  matrix  representation  for  (58)  is  given  by  (86). 

D.  Spheric!  Wave  Illumination 

The  final  vector  wedge  diffraction  problem  treated  here  is  the 
diffraction  of  a  spherical  electromagnetic  wave  incident  at  QE  along 

A  ^ 

the  direction  I  with  an  arbitrary  field  polarization  transverse  to  I. 
Such  a  spherical  wave  may  be  created  by  superimposing  the  fields  of 
z-directed  electric  and  magnetic  dipoles  at  s‘,  see  Fig.  5.  The  idea 
being  that  a  general  field  may  be  constructed  through  a  superposition 
of  the  TE  (to  z)  and  TM  (to  z)  waves  with  the  vector  potentials  T 
and  A,  respectively[19].  The  same  principle  was  used  in  the  pre¬ 
vious  problem  of  conical  wave  diffraction,  where  a  superposition  of 
electric  and  magnetic  traveling-wave  line  source  fields  generated 
an  arbitrarily-polarized  conical  wave. 

As  pointed  out  in  Chapter  II,  Section  C,  A  and  7  are  given  by 
m  =  zypegs(s,s')  and  F(s)  =  zepmgh(s ,  respectively. 

The  electric  dipole  moment,  zp05(  j7-s''|)  and  the  magnetic  dipole 
moment,  zpm5(!s-s'J)  generate  the  fields 
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*1  __  1  __ 

ne=  jj* V  x  ft  and  Em  =-  -V  x  F,  respectively. 

The  Green's  function  gs(s,?')  relates  the  field  at  i  due  to  a  source 

_  h 

at  s' , 

gs(s,s')  =  g(p,z;p* ,z';s“)  +  g(p,z;p' sz';e+), 
h 

where  g(p ,z;p' ,z' ;g+)  is  given  by  (62).  Now, 

Vs)  =  -v  x  zupegs(s,s')  =  -z  x  v  tpegs  (?,?') , 
and  similarly 

y?)  •  ;  x  vtpmgh(?,?'). 


It  can  be  verified  quite  easily  that 


(99)  -z  x  vtPe9s(s jS1 )  i  f 

,  i 


jkp  (s'-s  cos?) 


sin  eocot 


k(5l+s)C-j+j  s  s  s1niP-  (1  +  cos  5): 


(s'+s)‘ 


(s'+s)2  -  2s's  sin23Q  (1  +  cos  e) 


where  one  retains  only  terms  contributing  a  range  dependence  of 
0(l/p).  In  deriving  (99),  it  is  convenient  to  use  the  Green's 
function  in  (57)  rather  than  (62)  so  that  the  s/3p  and  the  1/p  3/3<j> 
operators  of  can  be  applied  directly.  Following  this  operation, 
the  transformation  p  =  s  sin  6Q  and  p'  =  s' sin  6Q  may  be  used 
along  with  the  approximation  of  (63);  then,  refining  the  proper 
order  of  the  range  dependent  terms,  one  obtains  (99). 
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Tha  complex  5  plane  topology  for  fU)  =  j(l  +  cos  5)  is  similar 
to  that  indicated  in  Fig.  3;  however*  the  integrand  of  ($9)  has  no 
branch  cut  singularities  in  the  complex  £  plane.  The  diffracted 
rays  corresponding  to  Re  may  be  obtained  via  the  saddle  point  ap¬ 
proximation  of  the  integral  of  (39).  Following  the  procedure 


outlined  in  Chapter  II,  the  diffracted  field,  may  be  expressed 

in  terms  of  I ( e+ ) ,  where 

jkp.  (s'-s  cosc)sinB.e"^s+s  ^ 


a  x  a  r 

<i>iu )  *-♦  j 


167r2jn  (s'+sJ^ss'siAjl+cosc) 
SDP(tt)  0 

sin2eo(l+cos5) 


cot 


«1) 
2n  J 


dr 


-  * 


r  jkpe  (s'-s  coss)sineoe"^s+s  ^ 
SDP(-tt)  16lf2  j  n  ( 1 s  *  ■ +S )  2-2ss ' : s i n2e0  ( 1+cos  5) 


cot 


(*£) 


sin2e  (1  +  cos  ?) 
e  5  +s  0  dc  . 


then 

(100)  -v  Jci(b')  -  l(g+)]  , 

c 

where 

e+  3  4  +  -j,1. 

V 

l(p+)  is  of  course  evaluated  via  the  modified  steepest  descent  method 
outlined  in  Appendix  1.  The  geometrical  optics  component  is  again 
obtained  from  the  residues  at  the  poles  of  the  integrand  of  (99). 

One  may  write  I(e+)  as  follows 
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(101)  *1(8)  *  +  (il<Pe|^r-  sin  60)^  cot(f$jFj 


rr~ 

s(s'+s) 


•  e'jks  FCKa+(e)]  - 


a /  e-jks'  \  , 

-VkpeTST  SineoJ^j 


-J 


•  3tr 


4njn  sin  er 


cot 


(s m 


S' 
s(s'+s) 


.  e"^ks  FLKa"(e)3  , 


where 


oo 


FCica-(e)3  =  2j 


^(g) 


•:Ka"(p)  f 


j  e"']T  dx,  and  k  =  ~p|  sin23c 


,  <a±(6) 


Rearranging  terms  in  (101),  and  utilizing  (100)  allows  one  to  write 

—  i  k  s  *  ■'  ■  » 

(102)  -  ♦  (jkpe  sin  S0)Dsj^5-  e-jks 


where  D$  is  given  by  (67). 


In  the  ray  optics  approximation,  one  may  use  El  =  -p 

6  ^  E 

A 

where  s  again  denotes  the  direction  of  the  diffracted  ray.  Using 


*  — A 

s  x  Hg 


S  X  <J>  =  -B  , 


(103) 


-  jks  * 


— <1  A  /  P  J1'3  ]  A  I 

£e  "  Bo(^““Pe  -51?'  s1n  Bo>  “sJstFnT  e'J 

-jks ' 


One  may  recognize  that  the  term  (ju)upQ  —x—rr-  sin  gj  represents  the 

c  ‘nfS  0 

A 

component  of  the  incident  electric  field  radiated  by  the  electric 
dipole  moment  evaluated  at  the  point  QE  on  the  edge. 


* j  .  Is"!  „  i  __  r !  j  il... 

rtQdlft,  C0  *  3q  ""  Q,,u  > 


(104) 


pd  A  r 

Ee  ”8o  E 


1  n  Is’ 

||  uSsjs (s 1 


+17 


-jks 
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From  duality,  one  may  obtain  -!*J  via  (100)  on  interchanging  H*j  by 

III  c 

Pe  Pm’  and  Ds  by  Dh’  respectively. 


(105) 


m 


<v 


-jks 1 

+^kpm  TiiT-  si"  6o>  Dhj 


S' 

s(s'+s) 


e-jfcs, 


where  is  identical  to  that  of  (67).  Alternatively,  (105)  may  be 
obtained  from  (|?xF),  where  it  is  understood  that  only  the  dif¬ 
fracted  field  contribution  of  (  ~  v  x  7)  leads  to  (105). 

Finally,  Hjj  =  (s  x  Ejjjx/eTu  Ip.  the  ray  optics  approximation, 
so  that 


_ jks  1  i  J  1  ■- ■  ■ 

(106)  -V  i0  (jUepm  ^  sin  90)  e-jks. 


-jks'  i 

One  may  recognize  ( jkp  sin  gQ)  in  (105)  to  be  E^  =  E^  •  <}>. 

-jks'  ^  ,,• 

Similarly,  (j(oepm  -y^-g-r  sin  &Q)  may  be  recognized  to  be  -H^  =-Hm  . 

A  -jks1  •  _•  a 

Also  (jkpe  sin  eQ)  in  (102)  is  equal  to  -H^j  =  -H^ 

Far-zone  conditions  for  both  the  source  and  the  field  point  are  as¬ 
sumed  in  the  above  analysis;  this  is  consistent  with  the  desired 
ray-optical  approximation.  One  may  re-write  (102),  (104),  (105) 


and  (106)  as 

r 


rd 

Ee  %  "8o  ^ 

■pd 


ri  n  s' 


«) 


-jks 


rd  :  ri  n  1  s1  -jks 

,  Era Ei  e 

f107)  1  nd  *  Hi  .  j  ?  -jks 

He  '  "♦  H»  Ds  j  iTs’+iT  e 

Tvd  *  ui  n  i  s7  -jks 

Hm  *  "eo  H-l  Dh  JsTs^+sT  e 
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The  total  diffracted  electric  and  magnetic  field  intensities 


are  respectively 

and 

e  m’ 

^  +  Hf . 

e  m 


In  terms  of  the  dyadic  diffraction  coefficient^ 


(108)  1 


-ShKt  e'JkS-  •*« 


Dr  =  -e'g  D 
E  o  os 


D 


In  the  convenient  matrix  notation,  (108)  is  given  by 


(109) 


rEdi 

K 

<v. 

o 

O  1 

in  _ 

Q  o 

1 

■ _  1 

'eJWe)' 

eJ(qe) 

k 

% 

'-Dh  0 
,°  -0, 

■nl(QE) 

,^(Qe) 

J 

1 


s' 

sTs+sTT 


s' 

s(s+'s‘) 


As  an  example  of  the  utility  of  the  above  results,  let  us  con¬ 
sider  the  diffracted  electric  field  of  an  arbitrarily-oriented 

A  ___  __ 

electric  dipole  moment,  u'  p  6(|s-s'|). 

Let  S'pe  =  Pep'o'  +  pe*'*'  +  pe2'*'>  where  V  V  and 

A  AAA 

Pez,  are  the  components  of  u!pg  in  the  p',<j>'  and  z‘  directions, 
respectively.  Table  I  gives  E,]  and  for  these  components. 
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The  diffracted  electric  field  is  readily  found  by  substituting 
E,],  E 1  into  (109). 

The  preceding  results  were  derived  for  0  <  <$>'  <  mr.  If  the  analyses 
are  repeated  for  <j>’  =  0,  nir,  i.e.,  for  grazing  incidence  along  the  sur¬ 
face  of  the  wedge,  in  each  case  the  diffraction  coefficients  Dh  are 
multiplied  by  a  factor  of  1/2  and  the  diffraction  coefficients  D$  =  0. 
This  also  may  be  deduced  by  considering  grazing  incidence  as  the  limit 
of  oblique  incidence.  At  grazing  incidence,  the  incident  and  reflected 
fields  merge,  so  that  1/2  of  the  total  field  propagating  along  the 
face  of  the  wedge  toward  the  edge  is  the  incident  field  and  the  other 
1/2  is  the  reflected  field. 


CHAPTER  IV 


DISCUSSION 


It  has  been  shown  that  when  electromagnetic  plane,  cylindrical, 
conical  and  spherical  waves  are  incident  on  a  perfectly-conducting 
wedge  the  resulting  diffraction  can  be  approximated  asymptotically 
in  the  form  of  the  geometrical  theory  of  diffraction.  Summarizing 
the  results  of  Chapters  II  and  III,  the  diffracted  electric  field 

(110)  E^s)  =  E1^)  •  FE(s,I)A(s)e”jks 


where  Ef1^)  is  the  incident  electric  field  at  the  point  of 
diffraction, 

A(s)  is  the  spatial  attenuation  which  describes  how  the  field 
intensity  vanes  along  the  diffracted  ray, 


(111) 


—  for  plane,  cylindrical,  and  conical  wave  incidence 

J~s 

(s  is  replaced  by  p  in  the  case  of  cylindrical  wave 
incidence) , 


(112) 


for  spherical  wave  incidence, 

\ 

De(s,I)  is  the  dyadic  diffraction  coefficient, 

%  =  D .(♦,♦' Mm'JbJ, 

L.  U  U  O  V  >•  v 
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(113)*  D  =  [{dV^FDcaV)]  *  €T(»“.n)FDca“(B")DJ 

s  sin  8q 

+{d+(s+1n)FCKa+{8+)]  + 


+  d“(8+,n)F[ica”(8+)]  }]  .where 


(114)  d±(8,n)  =  -  . . -  j  cot  ,  in  which  8  =  8+  =  (4>+<j>' ) 

nV2irk 

and 

(115)  FC»ca±(6)D  =  2j|^a*(8)l  ejKa"(e)  j  e"^2  dt. 

pu')| 

The  parameters  which  appear  in  F[tca±(8)3  are  defined  below 


(116) 


a+(8)  =  1  +  cos(-6  +  2nN+Tt) 


in  which  N±  is  the  positive  or  negative  integer  or  zero,  which  most 
nearly  satisfies  the  equations 


(117) 


2tlnN  -8  =  -ir  , 


2nirN  —6  =  it  . 


k  =  kL  is  the  largeness  parameter  in  the  asymptotic  evaluation  of 
the  pertinent  integrals  involved  in  the  formulation  of  the  dyadic 
diffraction  coefficient.  The  quantity  L  (appearing  in  k  =  kL)  may 
be  viewed  as  a  distance  parameter  which  depends  upon  the  type  of 
edge  illumination;  it  is  given  by 

s  sin280  for  plane  waves, 

(118)  L  =* jj+pT  for  cylindrical  waves, 
s's  si n^ 

— s  +~P~!L  ^or  conical  and  spherical  waves. 

s. 

-*For  grazing'  incidence  (^'=0,071)^ see  discussion  at  the  end  of  Ch,  III, 
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When  the  ray-fixed  coordinate  system  is  employed,  the  dyadic  dif¬ 
fraction  coefficient  may  be  expressed  as  the  sum  of  two  dyads.  One 
of  the  dyads  involves  D$,  the  scalar  diffraction  coefficient  for 
the  Dirichlet  boundary  condition,  and  the  other,  D^,  the  scalar 
diffraction  coefficient  for  the  Neumann  boundary  condition.  In 
turn  Ds,  depend  upon  trigonometric  functions  which  appear  in 

J.  .f. 

d^(g,n)  and  FCtca“(p)3»  which  involves  a  Fresnel  integral.  The 
latter  may  be  regarded  as  a  correction  factor  to  be  used  in  the 
transition  regions  of  the  shadow  and  reflection  boundaries. 

Outside  of  the  transition  regions  where,  tca'Ce)  =  kLa~(p)  >  10, 

F  is  approximately  equal  to  one.  Even  within  a  transition  region, 
usually  only  one  of  the  four  correction  factors  in  (113)  is  signifi¬ 
cantly  different  from  unity.  Curves  of  the  magnitude  and  phase  of 
F  as  a  function  of  kLa  are  presented  in  Fig.  6. 

Two  important  approximations  have  been  made  in  deriving  (110). 

In  accordance  with  the  geometrical  theory  of  diffraction,  it  has  been 
assumed  that  the  high-frequency  diffracted  field  propagates  along  it's 
ray  path  in  the  same  manner  as  the  geometrical  optics  field.  Con¬ 
sequently,  the  diffracted  electric  and  magnetic  fields  in  our 
asymptotic  solution  are  perpendicular  to  their  direction  of  propa¬ 
gation.  This  approximation  does  not  introduce  serious  error  if  both 
the  source  point  and  the  field  point  are  far  from  the  edge.  Secondly, 
in  deriving  (110)  it  has  been  assumed  that  kL  is  large.  However, 
based  on  the  extensive  numerical  study  of  asymptotic  solutions  of 
this  type  presented  in  Reference  9,  it  would  appear  that  generally 
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speaking  this  approximation  introduces  serious  error  only  when  kL  <  1. 
It  should  be  noted  that  the  accuracy  of  (110)  depends  on  both  approxi¬ 
mations,  and  a  study  of  the  error  introduced  by  the  two  approximations, 
when  either  the  source  point  or  the  field  point  (or  both  points)  are 
close  to  the  edge  has  not  been  carried  out. 

The  preceding  discussion  has  been  restricted  to  the  diffraction 
by  wedges  with  straight  edges;  however  the  geometrical  theory  of 

1  2 

diffraction  may  be  used  to  treat  the  diffraction  from  curved  edges,  * 
The  diffracted  ray  paths  are  determined  by  the  generalized  Fermat's 
principle  for  edge  diffraction,  and  the  conservation  of  power  flow 
in  the  resulting  astigmatic  bundle  of  rays,  see  Fig.  7a,  leads  to 
the  general  spatial  attenuation  factor 


(119) 


A(s) 


4 


s(Pc  +  s) 


where  the  caustic  distance  p  shown  in  Fig.  7a  is  given  by 


(120) 


1_ 

Pr 


S’ 


n  •  (I  -  s) 
pe  si"2|!o 


21 


for  spherical  wave  illumination  of  the  edge.  In  the  equation  above 
Pe  is  the  radius  of  curvature  of  the  edge, 

A 

n  is  the  unit  vector  normal  to  the  edge, 

A  A 

I,s  are  unit  vectors  in  the  directions  of  incidence  and  dif¬ 
fraction,  respectively. 


A  A 

60  is  the  angle  between  I  and  e,  the  tangent  to  the  edge  at 
the  point  of  diffraction  -  see  Fig.  7b. 

The  expression  for  A(s)  given  in  (119)  reduces  to  that  given  in  (111), 
if  p  is  set  equal  to  infinity  for  the  straight  edge;  furthermore, 

G 
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for  plane,  cylindrical  and  conical  wave  illumination  of  the  straight 
edge,  pc  =  ». 

The  diffraction  coefficient  is  assumed  to  be  independent  of 
Pe  to  a  first  approximation:  this  is  a  reasonable  assumption  because 
it  is  independent  of  the  curvature  of  the  incident  wavefront  to  this 
approximation,  as  has  been  demonstrated  here  in  Chapters  II  and  III. 
Furthermore,  the  validity  of  this  assumption  has  been  confirmed  in 
numerous  applications  of  the  geometrical  theory  of  diffraction  to 
structures  with  curved  edges.  Thus  in  accordance  with  the  postulates 
of  the  geometrical  theory  of  diffraction,  (110)  becomes 

(i2i)  r*(s)  =  r'(qE)  •  fE(s,i) J A  t  s)-  e-jks, 


where  is  given  by  Eq.  (112).  Alternatively,  if  the  incident  and 
diffracted  electric  fields  are  resolved  into  components  parallel 
and  perpendicular  to  the  planes  of  incidence  and  diffraction, 
respectively,  one  may  write  (121)  in  terms  of  matrix  notation  as 


0 


n 

+  s) 


By  introducing  the  proper  ray-fixed  coordinate  system  the 
polarization  effects  of  high-frequency  scattering  may  be  greatly 
simplified,  whether  this  involves  the  reflection  from  a  smooth 
curved  surface,  the  diffraction  from  an  edge,  or  the  diffraction 
from  a  smooth  curved  surface.  Specifically,  the  polarization  of 
the  scattered  field  may  be  related  to  the  polarization  of  the  in¬ 
cident  field  by  a  2  x  2  diagonal  matrix.  In  the  case  of  reflection 
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one  resolves  the  incident  and  reflected  fields  into  components 
parallel  and  perpendicular  to  the  plane  of  incidence,  as  is  well 
known,  and  in  the  case  of  diffraction  by  a  smooth  curved  surface, 
one  resolves  the  incident  and  diffracted  fields  into  components 
parallel  and  perpendicular  to  the  planes  of  incidence  and  dif¬ 
fraction,  respectively,  where  the  plane  of  diffraction  in  this 

case  contains  the  normal  to  the  surface  at  the  point  of  diffraction 

2 

and  the  diffracted  ray. 

In  summary,  employing  the  modified  Pauli -Cl emmow  method  of 
steepest  descent  and  introducing  a  ray-fixed  coordinate  system, 
a  compact,  dyadic  diffraction  coefficient  has  been  found  for  the 
perfectly-conducting  wedge.  This  diffraction  coefficient  is  valid 
in  the  transition  regions  of  the  shadow  and  reflection  boundaries 
for  a  variety  of  edge  illuminations.  From  a  practical  viewpoint, 
this  diffraction  coefficient  is  of  value  in  applying  the  geometrical 
theory  of  diffraction  to  antenna  and  scattering  problems  involving 
three-dimensional  structures  with  edges. 
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APPENDIX  I 

MODIFIED  STEEPEST  DESCENT  METHOD 
FOR  A  POLE  SINGULARITY  NEAR  THE  SADDLE  POINT 


Given 

(A-l)  I(k)  =  |  FjU)  eKfU)  <15, 
c 

and  the  contour  c  is  a  path  of  integration  in  the  complex  £  plane  such 
that  the  integral  converges.  The  following  procedure  outlines  the  ap¬ 
proximate  evaluation  of  1(0  for  large  k,  where  «  is  real. 

Let  f(c)  possess  a  saddle  point  at  5  =  £  (i.e.,  df/dc  =  0  at 
5  =  5S),  and  further,let  F^c)  possess  a  simple  pole  type  singularity 
at  c  -  £p.  The  procedure  to  be  outlined  allows  for  5  to  lie  in  the 
vicinity  of  £s. 

The  contour  c  is  deformed  so  that  it  passes  through  £s  along  the 
steepest  descent  path  (SDP).  If  any  singularities  of  the  integrand 
are  traversed  in  the  deformation  of  c  to  the  SDP,  their  contributions 
must  be  properly  accounted  for  in  the  evaluation  of  the  integral. 

The  SDP  in  the  complex  5  plane  is  given  by  the  equation 
(A-2)  Im  fU)  =  Im  fU$),  (subject  to  Re  f(0  <  0  for  |c|  °°). 


The  SDP  in  the  £  plane  may  be  mapped  into  a  straight  line  along  the 
Re  m  axis  in  the  complex  y  plane  via  the  transformation 


(A-3)  fU)  =  fUs)  -p  . 


Further,  the  descending  part  of  the  SDP  is  mapped  onto  the  positive  p 
axis  (this  determines  the  choice  of  the  branch  of  y  in  Eq.  (A-3).  Thus 


I(k)  =  e 


KfUJ 


F^U)  dy,  -A  <_  u  ±  B. 


-A 
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Assuming  that  is  analytic  in  the  nbhd  of  c$,  let 
(A-4)  F3(u)  =  Fj(0  ^  Cf(5)  -  f(«  )]  . 


F2(u)  is  now  analytic  in  the  nbhd  of  e  and  e  .  Thus, 

P  ® 


(A-5)  F2{p)  =  cmym. 


m-0 


and  is  a  valid  representation  for  F2(p)  about  the  saddle  point.  Let 


(A-6)  f(e)  -  fUp)  =  f(cs)  -  f(cp)  -  p2  =  -(u2  +  ja;, 


where 


(A-7)  a  =  j  [f(cs)  -  fUp)I]  . 


Using  (A-4),  (A-5),  and  (A-7)  in  I(<)  yields 

2 

»  r-  /  \ 

f 


,  ,,  <f<«s)  rB  Fa (ii)e~Kp  rfU  ) 

(A-8)  I  (k)  =  -e  — k -  dpw  -e 

-A  "  +ja 


.  c  |B£^e 

°i *2 ♦ *  ’ 


if  one  retains  only  the  first  term  in  the  MacLauren  expansion  for  F2(p). 
This  leads  to  the  first  order  saddle  point  result.  In  the  asymptotic 
approximation  one  may  extend  the  limits  of  integration  from  -«  to  «  y 
since  the  dominant  contribution  to  the  integral  is  from  the  immediate 
nbhd  of  the  saddle  point  at  p  =  0  when  <  is  sufficiently  large,  and 
elsewhere  the  integrand  falls  off  extremely  rapidly,  contributing 


negligibly  to  the  integral.  Therefore,  (A-8)  becomes: 

2 

r 


(A-9)  I(k)  *  -e  s  Cn\  §— - dp  . 


o  j  2 
—  00 


L  V  +  ja 
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FOa]  may  be  viewed  as  a  factor  which  accounts  for  the  effect  of  the 
pole  of  Fj(g)  which  may  lie  in  the  vicinity  of  the  saddle  point.  When 
!ca  co,  F[>a]-*»  1,  and  for  <a  >  10,  F[ica]  «1.  Thus,  when  the  pole 
singularity  is  sufficiently  far  removed  from  the  saddle  point,  FOa] 
may  be  replaced  by  unity,  and  (A- 10)  becomes 


(A- 12)  I(k)  %  e  s  Fj(5s) 

which  is  the  well  known  firs>t  order  asymptotic  result  for  the  case  of 
pole  not  close  to  the  saddle  point.  The  techniques  developed  in  this 
appendix  may  be  referred  to  as  the  "Pauli -Cl emmow  modification  of  the 
method  of  steepest  descent". 


76 


For  those  who  are  interested,  the  mapping  of  the  z  plane  onto 
the  p  plane  is  shown  in  Fig.  8.  Although  there  are  two  isolated,  first 
order  saddle  points  in  the  z  plane  (see  Fig.  3),  it  should  be  noted 
that  only  a  single  saddle  point  can  occur  in  the  p  plane.  Thus  only 
a  finite  region  of  the  z  plane,  adjacent  to  a  given  SDP,  is  mapped 
onto  the  p  plane.  In  the  present  case  there  are  two  SDP's  and  thus 
there  are  two  finite  regions  to  be  mapped  in  the  z  plane;  one  is 
defined  by  (0  £  u  £  2*,  -«  <  v  <  ~)  and  it  maps  onto  the  p  plane 
as  shown  in  Fig.  8a;  the  other,  defined  by  (-2tt  <  u  <  0;  -«  <  v  <  «) 
maps  onto  the  p  plane  as  shown  in  Fig.  8b.  In  carrying  out  the 
transformation  shown  in  Fig.  8a,  the  positive  branch  of  p  = 

J f(z)  -  f(g  is  used,  whereas  in  carrying  out  the  transformation 
shown  in  Fig.  8  b,  the  negative  branch  of  v~  Jf[z)  -  f  Us )  is  used. 
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[t  -  PLANE  REPRESENTATION  FOR  THE  REGION 

CORRESPONDING  TO  0£R«$£  2  it  (SEE  FIG.3) 
HO.  8a 


H  -  PLANE  REPRESENTATION  FOR  THE  REGION 

CORRESPONDING  TO  -  2  n  <  Re  |  <  0  (SEE  FIG.3) 
FIG.Sfe 


Fig.  8.  Complex  y  plane  topology. 
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APPENDIX  II 

FIELD  RELATIONS  FOR  SYSTEMS  WITH  NO  VARIATION  ALONG  THE 
z-  (AXIAL)  DIRECTION  IN  CYLINDRICAL  COORDINATES 


Let  v  =  vt  +  z  ~  t  where  vt  is  the  component  of  the  3-D 
Laplacian  operator  v,  which  operates  on  the  cylindrical  coordinates 

A 

transverse  to  z. 

Let  F  =  +  E"  and  IT  =  FF^  +  FFz,  where  F  and  FT  respectively 

are  electric  and  magnetic  field  intensities  satisfying  Maxwell  curl 
equations: 


(A— 13)  v  x  E  =  -joni  H  and  v  x  H  =  jue  E. 

E’t  and  Fft  are  the  transverse  (to  z)  field  components  and  E7,  Hz  are 
the  axial  fields  (along  z). 

A 

For  systems  uniform  in  z,  one  may  separate  out  the  z-dependence 

-jk  z 

in  E  and  H  by  assuming  an  e  dependence.  Next,  equating  the 
transverse  components  of  the  Maxwell  curl  equations  gives: 


(A— 14) 


vt  x  Ez  +  vz  *  Et  =  ”'JanJ  Ht 

VVVV  jue  Ft 


(where  v?  =  z  ~  ) . 

-jk_z  _ 

Finally,  using  the  e  L  dependence  in  E  and  FT  allows  one  to  solve 
for  FTt  via  (A- 14)  as  follows 

k2Ht  =  joie  Vt  X  Ez  +  Vz  X  vt  X  FTZ  +  Vz  X  v2  X  ff 
Or  p _  _  p  _ 

7t  X  Ez  -  jkz  vt  Hz  +  kz  Ht 
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i  .e., 


(A- 15) 


»t  Hz 


In  a  similar  fashion,  one  may  obtain: 

E  2  *  Vt  Hz 

(A— 16)  Et  =  -jkz  +  3co,i  p77T~ 


(Note  that  vt*  Ez  =  -z  *  vtEz  in  (A- 15).) 

Thus,  (A-15)  and  (A-16)  allow  one  to  calculate  Et  and  Ht  from 

a  knowledge  of  Ez  and  Hz  alone. 
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APPENDIX  III 

THE  AXIAL  FIELDS  OF  TRAVELING  WAVE  ELECTRIC  AND  MAGNETIC 
LINE  SOURCES  IN  THE  PRESENCE  OF  A  PERFECTLY  CONDUCTING 
INFINITE,  2-D  WEDGE* 


Case  1  (Traveling  wave  electric  current,  le  A  ) 

The  magnetic  vector  potential  A  associated  with  the  electric 
.  -jk  z' 

current  z  Ie  satisfies: 


(A-17 )  ft  =  z  Az,  and 

(vz+k2)Az  =  f  ie‘JkzZ'  dz.  _ 

«  P 

-jk  z 

Let  A  =  a  e  (since  the  source  and  the  field  A  must  have 
the  same  z-dependence  for  (A-17)  to  be  true).  Thus,  (A-17)  becomes 

(A-18)  (v?  +  (k^-k^)}a  =  -pi 

z  z  z  p 

Let 


(A-19)  k|  -  k2  -  k|  . 

The  fields  associated  with  Az  are  given  as: 

w  1  A 

(A-20)  H  =  magnetic  field  =  -  v  x  z  A  ,  and 
intensity  v 

v(v  •  A,z) 

(A-21)  E  =  electric  field  =  -jw  z  A,  +  — : - 

incensity 

The  boundary  conditions  on  Az  are  the  following: 
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2  g/\ 

(A-22)  E  |  n  n  =■  0  implies  ^  |-  A  n  n  =  0, 

p'<j>=0,niT  f  jwpc  3p  3z  4>=0 ,mi 


(A-23)  E  l  n  =  0  implies  joiA_  +  rr — 
z'4=0,nir  r  z  imiip 


Z  OQjpe  ^2  I^a0>nir 


For  (A-22)  and  (A-23)  to  be  simultaneously  true,  one  gets 

‘A-24>  AZ  W  *  0 

(and  therefore  az|^M„  -  0). 

In  addition,  Az  satisfies  the  radiation  condition  (for  a  time  dependence 
of  the  type  e+ju>t)  and  the  Mei xner  edge  condition.  It  follows  that 

(A-25)  az  =  ]il  Gg (p ,«f»;p *  ,<f. '  ;kt)  where 

(v2  +  k2)Gs(o,*;p'^'ikt.)  =  -  SW-'.M-tl). 


and  Gj  .  n  „  =0.  Also  satisfies  the  radiation  condition  and  the 
S'<f>=0,nTr  s 

Mei xner  edge  condition. 


Finally, 


i  ?  a 2  kf.  -jk  2 

E  =  — -  (r  +  )  A  =  -r^ —  a  e  1  H  -0  • 

z  June  9z2  z  jijpe  Z  ?  nz  * 


k;  -jk  z 
(A-26)  E_  =  le 

c  jwe 


Gs(ps*:p' j*1 ik^);  Hz  =  0. 
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,  “jk7z' 

Case  2  (Traveling  wave  magnetic  current,  Me  ). 

V 

The  vector  electric  potential  7  associated  with  the  magnetic 
*  -jk7z' 

current  zM  e  satisfies 
(A-27)  7  =  z  F  ,  and 

(V2  +  k2)F  =-e  f  Me"^2  dz> 

z  J  P 

••oo 

-jk  z 

Writing  Fz  =  f  e  as  before  * 

one  obtains: 

(A— 23)  (v2  +  k2)f  =  -eM  . 

i  t  z  p 

Thus, 

(A-29)  fz  =-eM  Gh(P,<j»;P,,V;kt) 


where 


(A-30)  (v2  +  k|)Gh(p,«j,;p*,4>';kt)  -  -  . 


The  boundary  conditions  on  Fz  (and  hence  fz)  and  are  identical 
and  are  obtained  from 


3H. 

_ _ i 

3<f) 


9H 


=  0  and  r~  =  0 

4>=0,nTT  3<?  l<j>=0,mr 


where 


H  =  -joj  F  +  ,  and  E  =  -  -  v  x  F. 

June  z 


It  follows  that  F  ,  f  and  satisfy: 
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3Gh  af  3F 

~  =  o,  r-r  =  0  and  ■  =  °- 

^  (j)=0,mr  ^  ij)=0,n7r  ^  <{)=0,nTr 

Finally  the  radiation  condition  and  the  Miexner  edge  conditions  must 

also  be  satisfied.  One  may  then  write 
l<2  ~jk  z 

(A-31)  Hz=  2  Gh(p’<,;p,»,J,,;kt)’  Ez  =  0* 
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